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Preface 


The  purpose  of  this  thesis  is  to  continue  investigation  of  the  applying  Porter's  high- 
gain  error  actuated  proportional  plus  integral  (PI)  design  techniques  to  the  development  of 
control  laws  for  an  Advanced  Tactical  Fighter  (ATF)  type  aircraft,  the  CRCA.  Initial  work 
was  accomplished  by  Capt.  Daryl  Hammond  and  the  results  of  applying  PI  control  based 
upon  output  feedback  were  very  promising.  The  results  obtained  of  this  effort  display  very 
robust  tracking  and  decoupling  of  outputs  with  the  fixed-gain  designs  and  a  much  more 
realistic  system  bandwidth.  More  investigation  of  the  adaptive  algorithm  will  be  necessary 
to  conclude  this  study  of  PI  control  for  the  set  of  inputs  and  outputs  selected  on  the  CRC.-\. 

.My  personal  thanks  and  gratitude  are  extended  to  all  who  have  helped  make  this 
thesis  effort  possible,  especially  my  thesis  advisor  Dr.  John  J.  D’Azzo  and  predecessor 
Capt,  Daryl  Hammond.  Thank  you  both  for  the  unending  support  ar.d  guidance.  The 
exceptional  support  and  resources  made  available  by  the  Flight  Dynamics  Laboratory  were 
also  greatly  appreciated. 

Finally,  I  would  like  to  thank  my  family  and  friends  for  their  kind  support  and 
understanding  through  this  long  and  difficult  endeavor. 


Jamie  Lynn  Foelker 
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Abstract 

Multivariable  control  laws  developed  by  Dr.  Brian  Porter  of  the  University  of 
''iilford.  England  are  used  to  successfully  perform  maneuvering  tracking  tasks  with  the 
NAS  A/Grumman  Control  Reconfigurable  Combat  Aircraft  (CRCA).  Porter's  method  is 
used  to  design  discrete  Proportional  plus  Integral  (PI)  control  laws.  Output  and  selected 
state  rate  feedback  are  used.  The  results  in  three  no  failure  flight  conditions  show  robust 
tracking  control  of  the  CRCA  for  five  selected  maneuvers.  Single  failures  are  introduced 
to  test  the  ability  of  the  fi.xed-g-'in  designs  to  successfully  control  the  aircraft  and  perform 
the  maneuvers.  The  time  responses  show  that  discrete  PI  control  law  can  make  the  CRC.A 
successfully  perform  all  five  maneuvers  for  two  of  the  three  control  surface  failures  investi¬ 
gated  in  two  of  the  three  point  designs.  The  step  response  PI  control  law  results  in  stable 
control  for  only  one  of  three  failu.-e  situations.  For  high  gains,  the  system  transfer  function 
becomes  asymptotically  diagonal  (the  outputs  are  decoupled).  Based  on  this  property,  the 
frequency  anaiysis  is  obtained  for  the  discrete  PI  design  using  each  output  with  respect  to 
its  associated  input.  Phase  margins  in  excess  of  ’5°,  gain  margins  of  greater  than  6dB,  and 
bandwudths  in  the  range  of  5-10  rad/sec  are  the  result.  The  adaptive  controller  displays 
a  larger  than  expected  roll  angle  output  in  two  of  the  maneuvers  as  compared  to  the  step 
response  PI  results.  An  adaptive  algorithm  using  a  recursive  least  squares  estimation  is 
run  with  failure  introduction  occurring  at  one  of  two  times  in  the  simulation.  The  adaptive 
rf'sults  also  display  decoupling  of  the  outputs  in  the  steady  state. 
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DISCRETE  PROPORTIONAL  PLUS  INTEGRAL  (PI) 
MULTIVARIABLE  CONTROL  LAWS  FOR  THE  CONTROL 
RECONFIGURABLE  COMBAT  AIRCRAFT  (CRCA) 

/,  [\TRODUCTlON 

1.1  DACKGROUSD 

Modern  aircraft  designs  present  the  control  engineer  with  great  challenges.  The 
demands  for  increased  maneuverability  can  be  met  only  at  the  expense  of  aerodynamic 
.^rability.  As  aircraft  become  highly  maneuverable  there  are  some  mission  fliglit  condition.s 
wlnrh  a  human  alone  cannot  control,  and  an  on-board  computer  must  be  used  for  basic 
.'Tability,  These  control  systems  are  termed  fly-by-wire  since  there  are  no  mechanical 
(■(mnections  between  the  pilot  and  control  surfaces;  only  electrical  signals  from  the  flight 
control  computer.  In  addition  to  the  role  of  stabilizing  the  aircraft,  the  flight  control  system 
must  also  provide  the  plane  with  the  ability  to  make  the  outputs  follow  input  commands. 
The  control  laws  developed  are  the  mathematical  equations  which  compute  the  required 
inputs  to  the  aircraft  or  plant. 

Modern  techniques  using  multivariable  control  have  a  distinct  advantage  over  their 
con'. muional  counterparts  in  that  only  one  control  system  need  be  developed  for  the  aircraft 
to  l)e  abie  to  perform  specific  tracking  tasks  where  more  than  one  set  of  inputs  and  outputs 
are  commanded.  Separate  loop  designs  are  necessary  in  order  to  make  each  output  track 
a  specific  input  when  single  input  single  output  (SISO)  designs  are  used.  This  can  make 
the  design  of  a  multiple  input  multiple  output  (MIMO)  system  extremely  involved  and 
tedious. 

.Another  issue  which  must  be  addressed  is  that  of  wdiether  the  designeil  lixed-gaiii 
flight  control  system  is  robust  enough  to  provide  adequate  control  when  the  aircraft  suffers 
a  control  surface  loss  or  malfunction.  This  problem  is  one  of  substantial  concern  to  the 
military.  The  first  most  impoitant  area  of  concern  is  that  of  returning  the  flight  crews 
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and  aircraft  to  safety  after  combat  damage  has  been  sustained.  Auoiiier  concern  is  iliat 
of  accomplishing  the  mission.  The  aircraft  must  remain  stable  and  retain  at  least  limited 
mission  capabilities. 

1.2  PROBLEM  STATEMENT 

Design  techniques  developed  by  Professor  Brian  Porter  of  the  L’ui\'crsity  of  Salfoid. 
f..ngland  allow  for  the  development  of  a  single,  fi.xed  gain  control  algorithm  which  U'cs 
output  feedback  to  make  the  plant  of  interest  perform  tracking  tasks  with  little  to  no 
coupling  of  the  outputs. 

The  goal  of  this  thesis  is  to  use  Porter's  multivariable  design  techniques  to  a.ssess 
the  effects  of  single  failure  conditions  on  control  algorithms  designed  for  a  healthy  aircraft 
performing  five  specified  maneuvers  in  three  different  points  on  its  flight  cmvelope.  The 
most  desirable  result  of  this  thesis  would  be  that  one  set  of  design  parameters  for  the  fi.Kcd- 
gain  designs  would  yield  adequate  stability  and  control  of  the  aircraft  in  all  three  flight 
conditions  as  well  as  provide  satisfactory  performance  when  a  single  failure  condition  is 
introduced.  Responses  obtained  from  single  failure  modes,  however,  prove  that  a  fi.xed  set 
of  design  parameters  is  an  optimistic  goal.  Thus,  the  introduction  of  adaptive  control  is 
deemed  necessary. 

This  thesis  effort  is  a  follow-on  to  the  work  accomplished  by  Capt.  Daryl  Hamnrond 
in  MULTIVARIABLE  CONTROL  LAW  DESIGN  EOR  THE  CONTROL  RECONFIG- 
URABLE  COMBAT  AIRCRAFT  (CRCA),  1988  (4).  Hammond’s  results  were  the  first 
obtained  using  Porter’s  control  techniques  on  the  CRCA  and  served  as  a  baseline  for 
performance  comparisons.  The  number  of  maneuvers  performed  will  be  increased  from 
Hammond’s  two  to  five  and  an  additional  failure  condition  will  be  simulated. 

l.-l  SUMMARY  OF  CURRENT  KNOWLEDGE 

Hammond  used  two  of  Porter's  design  methods  to  develop  Proirortional  plus  Integral 
(PI)  and  Proportional  plus  Integral  plus  Derivative  (PID)  control  algorithms.  Ilaininoiul 
used  four  different  techniques  to  design  PI  controllers; 
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1.  continuous  controller  based  on  known  plant  parameters. 

2.  discrete  controller  based  on  known  plant  parameters. 

3.  discrete  step  response  controller  b?«ed  on  unknown  plant  parameters,  and 

4.  adaptive  version  of  step  response  controller  to  take  plant  variations  into  considorat  ion. 

His  research  concluded  with  an  investigation  of  the  discrete  PID  ccmtiolim-  both  wiih  ami 
without  the  actuator  dynamics. 

Hammond  found  that  the  fi.xed.  high  gain  PI  confroller  ha.sed  on  t nowlcnige  of  tlw 
aircraft  parameters  displayed  robust  performance  in  normal  aircraft  flight  conditions  and 
that  a  fi.xed  set  of  design  parameters  provided  satisfactory  performance  in  70yc  of  the 
flight  configurations  tested.  He  also  found  that  the  discrete  PI  controller  outperformed 
its  continuous  counterpart.  However,  the  discrete  step  response  PI  controller  could  not 
adequately  control  the  left  flaperon  failure,  so  adaptive  control  was  applied.  The  adaptive 
controller’s  performance  was  e.xceptional  in  the  simulation  tested.  The  PID  controlh'r 
performed  very  well  in  simulations  where  the  CRCA’s  actuator  dynamics  were  not  included. 

1.4  .ASSU.MPTIONS 

The  following  assumptions  made  in  Hammond’s  thesis  will  be  \ised  in  this  effort,  and 
their  effects  on  the  thesis  scope  will  be  noted. 

1.  The  linearized  equations  of  motion  are  used  in  a  point  design  (initial  velocit}'  and 
altitude  are  assumed  fiAod  at  a  tiim  condition)  and  only  small  variations  arc  allowed 
about  the  trim  point. 

2.  The  outputs  used  are  assumed  to  be  measurable  and  readily  accessible. 

3.  The  aircraft  is  treated  as  a  rigid  body  so  that  bending  modes  do  not  appear  in  tlic 
linearized  equations  of  motion. 

4.  Fuel  consumption  is  considered  minimal  over  the  short  simulation  duratiotis.  .■^o  ma.s.s 
is  treated  as  a  constant. 
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5.  The  earth's  surface  is  used  as  the  inertial  reference  frame  which  can  again  be  at¬ 
tributed  to  short  simulations  and  aircraft  sensors  that  are  not  sensitive  enough  to 
detect  the  earth’s  rotation  (4:1-6). 

6.  The  atmosphere  is  assumed  fixed  to  the  earth’s  surface  and  wind  effects  are  not 
considered. 

/..5  APPROACH 

This  thesis  effort  focuses  on  the  two  methods  of  designing  discrete  non-adaptive 
PI  controllers  and  also  investigates  an  adaptive  approach  where  applicable.  Hammond's 
previous  work  is  used  as  a  baseline  and  the  new  research  has  been  built  from  his  existing 
computer  programs.  These  computer  programs  are  expanded  to  perform  the  additional 
flight  operations. 

The  CRC.-\  is  chosen  as  the  research  vehicle  primarily  because  there  is  a  groat  deal 
of  aerodynamic  data  available  for  single  and  multiple  control  surface  losses  and  failures. 
This  data  is  an  invaluable  aid  in  assessing  the  robustness  of  the  individual  PI  designs. 

Hammond's  design  parameters  are  used  initially  and  then  necessary  changes  made 
to  them  to  satisfactorily  control  all  outputs.  An  additional  goal  of  trying  to  obtain  “rea¬ 
sonable"  controller  gains  is  added.  Frequency  analysis  will  be  performed  to  ensure  that 
the  phase  and  gain  margins  specified  by  MIL-F-9490C  are  met  and  a  realistic  bandwidth 
is  obtained. 

The  research  will  be  performed  on  the  powerful  software  package  MATRIX  y-  .\H 
necessary  modifications  have  been  made  to  the  computer  programs  in  order  to  run  the 
siiiiulations  on  Version  7.0.  The  System  Build  feature  of  MATRIX^  (6)  allows  for  straight¬ 
forward  multivariable  design  and  analysis. 

1.6  OVERVIEW 

The  remainder  of  this  thesis  is  organized  according  to  the  following  schedule.  The 
aircraft  models,  actuators,  and  ARMA  representation  are  discussed  in  Chapter  2.  The 
theory  behind  the  different  control  law  design  methods  are  presented  in  Chapter  3.  The 


1-4 


The  results  obtained  from  tlie  discrete  PI  controller  can  be  found  in  (Chapter  1,  aiui  the 
step  response  and  adaptive  results  are  contained  in  Clia[)ter  5.  (.'hapter  o  contains  a 
summary  of  the  results  obtained  and  includes  recommendations  for  future  work  u^ing 
Porter's  methods  on  the  CRCA  flight  control  system. 
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II.  DESCRIPTIOS  OF  MODELS 


J.  1  .1  IRCR.AFT  DESCRIPTIOS 

The  Control  Recoafigurable  Combat  Aircraft  (CRCA)  is  a  NASA/Grumnian  Ad¬ 
vanced  Tactical  Fighter  ( ATF)  type  design.  The  CRCA  was  designed  for  high  performance 
with  relaxed  static  instability  of  12%  at  low  airspeeds  (12:3-1).  The  plane  is  not  slated 
for  production,  but  it  is  being  used  to  assess  combat  damage  survivability.  The  CRC.A. 
shown  in  Figure  2.1,  has  nine  separate  control  surfaces  consisting  of  two  canards,  a  i  iiddvr. 
and  two  flaperons  and  an  elevator  on  the  trailing  edge  of  each  wing.  The  three  siirfaci's 
on  the  trailing  edge  of  each  wing  are  commanded  together  in  this  thesis  effort  in  ‘Tdor 
to  preserve  the  laminar  airflow  over  the  wings.  The  canards  have  30°  of  positive  dihedral 
(upward  cant  from  the  horizontal)  and  are  '•all-flving”,  meaning  that  the  longitudinal  and 
lateral-directional  equations  of  motion  for  the  aircraft  are  coupled. 

2.2  FLIGHT  PHASES 

The  aerodynamic  data  on  the  CRC.-\  is  available  for  four  design  points  and  the  four 
design  points,  or  flight  conditions  are  representative  of  the  aircraft’s  performance  range. 
Descriptions  of  the  four  flight  pha.ses  w'as  obtained  from  .Appendix  C  of  Reference  ( 12)  and 
Hammond’s  Chapter  2. 

The  first  flight  condition  is  Air  Combat  Maneuvering  (ACM)  Entry.  This 
flight  phase  is  representative  of  aircraft  conditions  normally  encountered  at  the  initiation 
of  an  air-to-air  engagement  or  weapons  delivery.  .Altitude  can  vary  from  10,000  to  .30.000 
feet  with  velocities  just  under  Mach  1. 

The  next  flight  phaise  is  Terrain-Following/Terrain-Avoidance  (TFTA).  This  is 
a  low  altitude,  high  velocity  flight  condition  commonly  used  to  avoid  enemy  radar  detection. 
Velocities  are  as  high  as  those  used  in  .ACM  Entry,  but  with  altitudes  as  low  as  200  feet 
above  local  ground  level.  The  dynamic  pressure,  7,  is  very  high  for  this  condition  and 
rapid  changes  in  pitch  angle  can  be  expected. 

The  third  point  design  is  ACM  Exit.  This  condition  is  very  difficult  to  control  since 
the  aircraft  has  low  airspeed  with  a  high  angle  of  attack  (~  30°),  steep  bank  angle  (~  70°  !, 
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and  a  high  load  factor  (~  2g).  The  aircraft  has  been  pushed  very  close  to  its  limits  by  the 
end  of  the  air-to-air  encounter. 

The  fourth  flight  phase  is  Short  Take-OfT  and  Landing  (STOL).  This  condition 
has  the  lowest  airspeed  of  all  four  and  is  the  most  critical  flight  phase.  The  lauding 
requirements  are  a  ma.ximum  landing  distance  of  1-500  feet  and  a  crosswind  of  .'50  knot.s. 
This  data  was  unfortunately  not  available  in  time  to  be  incorporated  in  this  thesis  effoit. 

Table  2.1  (4:2-3)  summarizes  the  characteristics  of  the  C’RCW's  flight  phases.  Tlie 
point  designs  accomplished  in  this  thesis  are  valid  since  the  maneuvers  performed  do  not 
stray  significantly  from  the  design  points. 


Table  2.1.  Flight  Conditions 


.MISSION  SEG.MENT 

ALTITUDE.  FT 

1  MACH  NO. 

LOAD  FACTOR 

] 

ACM  ENTRY 

30,000  1 

0.9 

.  ’  1 

i 

TF/TA 

SEA  LEVEL 

0.9 

ACM  EXIT 

10,000 

0.275 

STOL 

1200 

0.185 

3.3  DETERMINATION  OF  AIRCRAFT  MODELS 

2.3.1  STATE  SPACE  MODELS  Linearized  state  space  models  of  the  CRC.A  have 
been  obtained  utilizing  computer  programs  at  the  Flight  Dynamics  Laboratory  at  Wright 
Patterson  Air  Force  Base  (4,  7).  .■Ml  aircraft  models  used  in  this  thesis  are  those  obtained 
from  Hammond,  with  the  exception  of  the  25%  rudder  failure  case  models  which  are 
extracted  from  Capt.  Kurt  Neumann’s  thesis  computer  files  (7). 

The  result  of  using  the  CRCA’s  nonlinear  equations  of  motion  and  control  and  sta¬ 
bility  derivatives  is  a  linearized  model  containing  nine  states  and  nine  inputs,  as  shown  in 
Table  2.2  (4:2-9)  for  the  .■\CM  Entry  flight  condition. 
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Table  2.2.  Aircraft  State  Space  Matrices  for  ACM  Entry 


.0119 

-.0186 

-31.2350 

-32.1804  .0000 

.0000 

.0000 

.0000 

.0000 

.0324 

-1.0634 

894.4548 

-1.0634  .0000 

.0000 

.0000 

.0000 

.0000 

.0002 

.0069 

-.6015 

COOO  .0000 

.0000 

.0000 

.1)000 

,0000 

.0000 

.0000 

1.0000 

.0000  .0000 

.0000 

.0000 

0000 

oooo 

.0000 

.0000 

.0000 

.0000  .0000 

.0000 

.0000 

1.0000 

,0000 

.0000 

.0000 

.0000 

.0000  .0000 

.0000 

.0000 

1.0000 

.0349 

.0000 

.0000 

.0000 

.0000  .0000 

.0360 

-.0929 

.0349 

-.0994 

.0000 

.0000 

.0000 

.0000  .0000 

.0000 

-27.8066 

-2.0376 

.4913 

.0000 

.0000 

.0000 

.0000  .0000 

.0000 

2.4582 

-.0241 

—  .4377 

.0411 

.0411 

.1322 

.0866  .1322 

.0866 

.1018 

.1018 

.0000  ■ 

-.3163 

-.3163 

-.9597 

-.6194  -.9597  - 

.6194 

-1.0183 

1 

o 

<o 

.0000 

.1014 

.1014 

-.0284 

-.0215  -.0284  - 

.0215 

-.OjoO 

-.0200 

.0000 

.0003 

-.0003 

-.0002 

—  .0001  .G0o2 

.0001 

-.0001 

.0001 

.0006 

.0000 

.0000 

.0000 

.0000  .0000 

.0000 

.0000 

.0000 

,0000 

.OOUU 

.0000 

.0000 

.0000  0000 

.0000 

,0000 

.otiuO 

.0000 

.0000 

.0000 

.0000 

.0000  .0000 

.0000 

,0000 

.01100 

.0010 

0762 

-.0762 

.2219 

.2011  -.2219  - 

.2011 

.1109 

-.1109 

.1141 

.0486 

-.0486 

.0029 

.0021  -.0029  - 

.0021 

.0021 

-.0021 

-.0.544 
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The  eight  states  As  mentioned  previously,  the  three  trailing  edge  control  surfaces  on 
each  wing  are  commanded  together.  This  effectively  reduces  the  number  of  control  surfaces 
from  nine  to  the  following  five; 

Sci  =  left  canard 

6cr  =  right  canard 

Stfi  =  left  trailing  edge  flaperon 

Ster  =  right  trailing  edge  flaperon 

Srud  =  rudder 

The  columns  of  the  corresponding  left  and  right  trailing  edge  .surfaces  of  the  B  matrix 
are  .summed  to  form  a  reduced  input  B  matrix.  Now  that  there  are  only  five  ini)uts.  only 
five  outputs  can  be  controlled  using  output  feedback.  The  five  chosen  outputs  are 

V  =  forward  velocity 
3  -  sideslip  angle 
=  pitch  angle 
<P  =  bank  angle 
r  =  yaw  rate 

The  9x9  A  matrix  does  not  have  full  rank  since  the  state  ip  is  redundant.  Therefore, 
the  fifth  row  and  column  are  removed  and  A  is  reduced  to  an  eight  state  model.  .Also, 
Hammond  found  that  the  remaining  eight  states  had  to  be  reordered  so  that  .Autoregressive 
Moving  .Average  (.ARM A)  models  could  be  generated  more  directly.  Table  2.5  (4:2-10) 
contains  the  resulting  ACM  Entry  model  of  the  CRCA  used  in  this  and  Hammond’s  thesis. 
The  order  of  inputs  and  states  remains  the  same  for  the  rest  of  the  flight  conditions  used 


where  the  eight  remaining  states  are 


u  =  velocity  in  aircraft  X  axis 
w  =  velocity  in  aircraft  Z  axis 
q  =  pitch  rate 
0  =  pitch  angle 

^  =  yaw  angle 

0  =  roll  angle 

J  =  sideslip  angle 
p  =  roll  rate 
r  =  yaw  rate 


The  aircraft  open  loop  eigenvalues  for  the  three  point  designs  used  are  listed  in 
Table  2.3  and  the  transmission  zeros  in  Table  2.4. 

Table  2.3.  Open  Loop  Eigenvalues  of  8x8  CRCA  Models 


ACMENTRY 

ACMEXIT 

TFTA 

-0.0078±j0.0495 

-0.0562 

1.6669 

-0.2482±jl.8323 

-2.0156 

-3.3281 

0.0220±70.0363 

-0.2139 

-0.3059 

0.4718 

-1.7162 

-0.6572±jl.7084 

0.0236 

-0.0457 

-0.0583 

-0.7559±;2.9853 

-2.4418±;4.4400 

-5.4284 

Table  2.4.  Open  Loop  Transmission  Zeros  of  8x8  CRCA  Mouels 


ACMENTRY 

ACMEXIT 

TFTA 

-1.6844 

-0.0902 

-2.8974 
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Table  2.5.  Equivalent  Aircraft  8x8  Matrices  For  ACM  Entry 


.0000 

.0000 

0000 

.0000 

1.0000 

.0000 

.0000 

.0000 

.0000 

.0000 

.0000 

.0000 

.0000 

.0000 

1.0000 

03.50 

-32.1804 

.0000 

-.0119 

-.0186 

-31.2350 

.0000 

oooo 

.0000 

-1.06.34 

.0000 

-  0324 

-1.0634 

894.4548 

.0000 

OdiiO 

0000 

.0000 

.0000 

0000 

.0069 

—  .6015 

,0000 

,01)00 

0000 

.0000 

.0360 

0000 

0000 

.0000 

-.0929 

11.349 

-  '.'OO-l 

.0000 

.0000 

0000 

.0000 

.0000 

—  27,8066 

-2  0.3:1) 

.0000 

.0000 

0000 

0000 

.0000 

2.4.582 

-  0241 

1377 

.0000 

1  .0000 

.0000 

.0000 

.0000 

.0000 

1  .0000 

.0000 

.0000 

.0000 

.0411 

.0411 

.3206 

.3206 

,0000 

-..3163 

-.3163 

-2.5974 

-2.5974 

.0000 

.1014 

.1014 

-.0699 

-.0699 

.0000 

,0003 

-.0003 

-.0004 

,0004 

.0006 

.0762 

:  -.0762 

.5339 

-.5339 

.114  1 

.0486 

-.0486 

.0071 

-.0071 

-.0544 

x=  [ 

0  0  u 

W  q  3 

P  r  ]' 

1  T 

u.  = 

^tel  ^ter 

^rud  j 

The  output  matrix 

is 

'001 

.0.349  0 

0  0  0  ' 

0  0  0 

0  0 

1  0  0 

'  = 

1  0  0 

0  0 

0  0  0 

0  1  0 

0  0 

0  0  0 

0  0  0 

0  0 

0  0  1 
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Simulations  will  be  performed  in  the  three  flight  phases  ACMENTRY,  ACMEXIT. 
and  TFTA  and  the  following  single  failure  cases  are  investigated  in  the  ACM  Flntry  and 
TFTA  flight  modes: 

1.  30%  loss  of  left  trailing  edge  (combination  of  three  surfaces) 

2.  '0%  loss  of  left  canard 

3.  2')%  loss  of  rudder 

No  failures  are  investigated  in  ACM  Exit  since  the  flight  condition  is  very  difTicult  to 
control  with  the  healthy  plant.  .Appendix  A  contains  the  reduced  order(SxS)  state-space 
models  for  all  flight  conditions  and  for  the  single  failures  used  in  this  thesis. 

2..C2  ACTOREGRESSlVt:  MOVISC.  AVERAGE  (ARMA)  MODEL  dokor  and 
Keviczky's  (  1)  technitpie  is  used  to  transform  a  plant  model  from  state  space  fortn 

X  =  .lx  +  Ru 
y  -  Gx 


to  an  .ARM.X  representation  where 

^  A^yik  -  l)r  +  •■  ■  +  A.syik  -  y  )T  =  B^u{k  -  ])T  +  ■  ■  ■  +  Byu{k  -  S  )T  (2..3) 

and  the  matrices  .li  through  .4,v  and  B\  through  B\  are  the  .ARM.\  coefficients.  I  his 
equation  is  used  by  Porter  and  Fripp  (9.  8)  with  no  noise  source  considered.  The  reduced 
order  .AR.MA  model  is  determined  using  the  ratio  of  the  number  of  states  n  to  the  number 
of  inputs  m,  y  ~  n/m.  For  the  eight  state,  five  input  model  used  in  this  thesis,  two  .1, 
and  B,  .AR.M.A  coefficients  an  sufficient  to  form  the  output  Equation  2..') 

F.fjiiation  2. .3  ran  also  be  expressed  as 


(2.(;) 


y(kT)  =  e^(kT)o{kT)  =  f(kT)o^[kT) 


where,  for  ease  of  simulation, 


’  .4i  ■ 

-y(k  -  1  )T 

-4,v 

1 

1 

^(A;T)  = 

u{k  -  1  )T 

.  . 

u[k-N)T 

Karh  coefficient  of  6{kT]  is  a  column  vector  containing  the  P  terms  of  that  ARMA  coeffi- 
cierit.  Thus,  the  elements  iii  the  rows  of  .4)  are  transposed  and  placed  in  one  column  that 
contains  the  values.  Then  the  elements  of  .42  occupy  the  i.:.xt  P  positions.  Similarly, 
the  elements  of  /fj  and  Bi  each  occupy  P  positions  in  6[kT). 

f.ach  term  of  o(kT)  is  an  \  m  block  diagonal  sub  matrix.  The  past  outputs 
iy(k  -  i)T,  with  1  <  t  <  .V,  for  each  of  the  m  inputs  appear  as  diagonal  blocks.  The  ini)uts 
n(k  -  t)T.  with  1  <  t  <  -V,  appear  in  diagonal  subblocks,  each  of  size  m  x  1,  with  only  ii, 
appearing  in  each  subblock.  See  the  example  matrices  on  the  next  page. 

The  reduced  order  model  for  the  CRC.4  uses 

OikTf  =  [aJ  aI  bJ  e  (2.7) 

Oiklf  =  [-.yT'{k-l]T  -y^(fc-2)T  1)T  u^(fc-2)r]^G  /?'’^‘“(2.8) 

The  program  used  to  generate  the  .\RMA  coefficients  for  each  flignt  condition  along  witli 
the  coefficients  for  each  flight  condition  is  listed  in  ,\ppendix  B.  The  .\RMA  representation 
of  the  system  is  used  to  run  the  adaptive  algorithm  of  references  (8,  1.3). 
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I'or  >'xa:;ri)l>\  witli  n  =  2  with  2  inputs  and  outputs: 


.dl(l.  1) 

-  1) 

0 

■di(1.2) 

-t/2(fc  -  1) 

0 

.Ai(Tl) 

0 

-t/il/-’  -  1) 

.d,CT2) 

0 

-  1/2  (  -  1 ) 

■Mil.l) 

-  2) 

0 

.-12(1.2) 

-y2ifc-2) 

0 

■■l2(2.  I ) 

0 

-U\(^  -  2) 

.12(2,2) 

0 

! 

I 

o  = 

i) 

ui{k  -  1 ) 

0 

Bi(1.2) 

ii2{k:  —  1 ) 

0 

/lii2.  1) 

0 

ui(k  ~  1) 

/ill 2.  2) 

0 

U2{k  ~  1 ) 

/i2(l.  1) 

uxik-2) 

0 

Ih{E2) 

U2{k  -  2) 

0 

Ih'M-A) 

0 

ui{k-2) 

/i2(2.2) 

0 

U2{k~2) 

J.J.J  .\<'"r!.'.\  'TOR  MODELS  The  CRCA  has  no  mechanicai  coimcctions  botwoeii 
rln'  pilot  and  tiif  control  surfaces  and  is  thus  termed  a  fly-by-vvire  control  .system.  In  ortier 
U)  a(id  vnlidity  to  the  designs  accomplished,  actuator  dynamics  must  be  taken  into  consiti- 
l■,";lIi<Jn.  In  addition  to  actuator  dynamics,  the  physical  limitations  of  maximum  delieciions 
and  d'dl.Ttion  rates  must  also  be  incorporated  into  the  actuator  mod<  ■1.  Table  2.fi  li.-^ts  the 
r;iiiitatif)n,..  for  the  three  diiferent  control  stirfaces. 


Ha.mmorifl  wtjrked  pritiiaril v  with  tin'  fir.st  order  actuator  mode] 


20 


I.S  +  20) 
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Table  2.6.  Control  Surface  Position  and  Rate  Limits 


Surface 

Position  Limit 
(deg) 

Rate  Limit 
(deg/sec) 

Canards 

+60  -30 

±  100 

Trailing  Edges 

+30  -30 

±  100 

Rudder 

+20  -20 

±  100 

second  order  approximation  of  the  actuator  model  is  used  exclusively  throughout  this 
thesis  effort. 


•^con(ro/('^) 


5113 

(s  +  53±  j48) 


(2.10) 


.A.  description  of  the  CRCA  is  described  by  a  set  of  linear  differential  equations 
and  an  eight  state  model.  The  aircraft  models  are  linearized  about  3  flight  conditions, 
TF/TA,  .ACM  Entry,  and  ACM  Exit  and  include  second  ord  'r  actuator  dynamics.  The 
autoregressive  moving  average  (ARM.A)  representation  of  aircraft  is  also  presented.  The 
PI  control  law  development  is  presented  in  Chapter  3. 
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III.  DESIGN  PROCEDURES 


3.1  OVERVIEW 

The  PI  controller  design  techniques  developed  by  Porter  are  implemented  using  out¬ 
put  feedback.  Using  output  feedback  is  often  more  advantageous  than  u.siiig  state  feedback 
if  a!i  of  the  states  arc  not  readily  available  and  estimates  would  have  to  be  used  ( l. 
Using  Porter's  method,  tracking  can  be  achieved  with  little  to  no  coupling  occurring  Im- 
tween  outputs.  The  use  of  high  gain  in  the  forward  loop  results  in  an  asymptotic  transfer 
function  for  the  closed  loop  system  composed  of  fast  and  slow  modes.  “The  slow  modes  are 
asymptotically  uncontrollable  or  unobservable,  and  thus  the  output  response  is  dominated 
by  the  fast  modes.”  (3:661) 

3.3  DISCRETE.  FKXED-GAIS  FI  CONTROLLER  THEORY 

The  plant  must  satisfy  four  necessary  criteria  before  any  of  Porter's  design  techniques 
can  be  used: 

1.  The  plant  must  be  completely  controllable. 

2.  The  plant  must  be  completely  observable. 

B  .\ 

i.  rank  —  n  +  I 

D  C 

T  Transmission  zeros  can  be  computed  from 

si  -  A  B 

=  0 

I  -C  D 

for  the  completely  controllable  and  completely  observable  plant.  .\’o  transmission 
zeros  can  be  located  at  the  origin,  or  the  plant  is  functionally  nnconli'ollable.  1  lie 
number  of  transmission  zero.s  for  the  plant  where  I  —  m  is  n  —  m  —  d,  where  il  is  the 
rank  deficiency  of  CB  (3:668). 


The  continuous  aircraft,  or  plant,  may  be  represented  by  equations  of  the  form 


Xi(t) 

•111  -'ll  2 

^l(0 

0 

= 

+ 

X2{t) 

.421  -■I22 

X2(t) 

B2 

C,  C2 


where  the  dimensions  are 


ii(/) 

X2{t) 


X\  :  ( n  —  m )  X  1 

X2  :  TO  X  1 

111  :  {n  -  m)  X  {n  —  m) 

.422  :  TO  X  1 

Bi  ■  rn  X  TO  and  has  rank  m 

C2  :  m  X  TO  and  has  rank  to 


(3.1) 


(3.2) 


Notice  that  the  B  matrix  in  Table  2.5  and  the  B  matrices  for  the  rest  of  the  CRC.A's 
flight  conditions  are  not  of  the  required  form  in  Equation  3.1.  A  transformation  could  be 
performed  to  express  the  CRCA  state  space  matrices  in  the  form  of  Equation  3.1.  However, 
this  step  has  been  proven  to  be  unnecessary  (3:668-669),  and  the  design  can  be  performed 
for  any  form  of  the  state  equations. 

The  rank  of  the  matrix  product  CB  (the  first  Markov  parameter)  must  be  checked. 
For  a  plant  that  has  rank(C5)=m,  the  plant  is  termed  “regular”  and  only  requires  output 
feedback.  For  a  plant  like  the  CRCA  where  rank(CB)<  m,  the  plant  has  a  rank  deficient 
first  Markov  parameter  and  rate  feedback  of  the  state  vector  i\  must  augment  the  output 
feedback.  Figure  3.1  illustrates  the  control  law  used  for  irregular  plant.s. 

The  feedback  signal  with  rate  feedback  incorporated  is  now 

w  =  y  +  Mil 
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Figure  3.1.  Discrete  PI  Controller  -  Irregular  Plant 

(4) 

From  Equations  3.1  and  3.2  it  can  be  shown  that 


\ 

w  = 

c, 

C2 

Xl 

■4  12 

1 

X2 

1 

12 

— 

F2 

•Tl 

^2 

From  Figure  3.1  it  can  be  seen  that 


e  =  r  —  w 


w  here  r  =  commanded  input  signal. 

E.xpressing  the  state  equation  in  the  discrete  time  domain  yields  the  difTcrencc 

tion 


x(/t  +  1)T  =  4>{T)x{kT)+  'li{r}u{kT) 


where 


^T)^ 


,AT  _  ^12{T) 

<t>2l{T)  4>22{T) 


fT 

^{T)  =  /  e-^^dr 
Jo 


MT) 


anu  ihe  ainerencc  equauuu  re^reseulliig  Lhe  luiegral  of  the  error  is 


;(it  +  l)T  =  zikT)  +Te{kT) 


=  z{kT)  +T[r{kT)  -  w(kT)] 

=  z[kT)  +Tr{kT)  -  TFxx^{kT)  -  TF2X2{kT) 


The  composite  difference  equation  for  the  open  loop  of  Figure  3.1  is 


z{k  1  )T 
Xi{k  +  1  }T 
X2(k  +  1  )T 


It  -TFi  -TF2  z(kT)  0 

0  4>n  4>i2  X\{kT)  +  tf'i  u[kT) 

0  021  022  ^2{kT)  IP2 


+  0  r(A:T) 


z(kT) 

y  =  [  0  C,  C2  ]  ii(fcT) 
^2ikT) 


where  the  control  law  is 


u{kT)  =  -[Kte(kT)  +  K2z(kT)] 

=  Yi^ur(hT)  -  A'.F,x,(fcT)  -  h\F2X2{kT)  +  K.zikF 
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The  closed  loop  equation  for  the  system  can  now  be  written  as 


c{k  -t  i)T 

tc 

F-h 

1 

bn 

t 

t _ 

z(kT) 

xi(k  +  l)T 

= 

d>u  —  7  ti’i  A 1  Fi  Oi2-7d’iAiF2 

Xi(kT) 

X2{k  +  l)T 

^!^’2A2  d>21  —  J>l’2T\Tl  022  -  f^’2l'^\F2 

X2{kT) 

Tit 

jd’i  Ai 

T  *^'  2  A  1 


rlhT) 


No  changes  appear  in  Equation  3.9,  the  system  output. 


(3.10] 


The  form  of  Equation  3.10  that  is  in  Porter's  paper  (2)  is  obtained  by  making  the 
substitutions 


<^(T]  =  =  I  +  AT  + 


.[2t2 

2! 


+  . .  . 


(3.11) 


Hm^e-''^  =  /  +  .4T  +  0(H.O.T.) 


lim  -I]  =  /1  +  0(H.O.T.) 

t^oT 


0\\ 

012 

A\\T  -t-  Iji-t 

A12T 

4>12 

022 

A21T 

A22T  4-  It 

0 

1 

0 

0 

fu'i(r) 

=  ;p{T) 

= 

~ 

Jo 

B2 

1 

B2 

B2 

fd-2(T) 
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Equation  3.10  can  now  be  written  as 


:(k  q-  1 ) T 

It  -TFi  -TF2 

x[kT) 

xi{k  -f-  l)T 

= 

0  •dllT'  +  Iri-t  Ax2T 

xiikT) 

X2(k  +  1)7" 

52A2  ■'^2\T  —  B2h\Fi  A22T  F  It  — 

X2[kT) 

Th 

0 

02  A  1 


r[kT) 


(3.13) 


The  output  equation  has  not  changed  from  Equation  3.9. 

Singular  perturbation  methods,  where  f=j  — *■  oo  are  the  perturbation  variable,  is 
used  to  block  diagonalize  the  closed  loop  system  of  Equation  3.13  and  to  determine  the 
fast  and  slow  modes  of  the  system  (3,  2). 

Equation  3.13  must  first  be  represented  by  the  form 

r(kT)  =  Ax{kT)  +  Br{kT)  (3.11) 


xi{k  -h  1  }T 

1 

x(kT)  + 

= 

X  2  ( ^'  "f"  1 )  r 

•da 

B2 

yi^T)  = 


Cl  C2 


x{kT) 


where 


XxikT) 


and  the  matrix  dimensions  are 


x{kT) 

xx(kT) 


■4i  :  (n  -  m)  X  {n  -  m) 

A4  :  m  X  m 

82  '-  m  X  m  and  lias  rank  m 


(3.1.5) 
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A  transformation  matrix,  T,  can  be  used  to  transform  Equation  3.14  into  bl(;ck 


diagonal  form  (3) 


i,(k  +  i)r 

A,  0 

xJkT) 

1 

+ 

r  B, 

Xf(k  +  1)T 

0  A/ 

_  j 

1  . 

r{kT) 


y{kT)  = 


C,  Cj 


x,(kT) 

XfikT) 


D’Azzo  explains  the  steps  to  block  diagonalize  Equation  3.14  (3;G63-G67). 
D'Azzo’s  results  to  those  that  Porter  obtained  for  the  discrete  plant  (2:1210), 


Aj  —  I  T Ao 


-A'f'A'2 

Ai2  r '  2  A 11  -  .d  1 2  A" 


.4/  =  /  +  .44 


A4  -  -  52/11^2 


B,  =  TBo 


Bo 


0 

Al2i^2'' 


B I  =  B2  A 1 

C,  =  Co 


Co 


C2f2"‘A'f’A'2  Ci-CCAf'/^ 


C/  =  C2 


(3.16) 


(3,17) 

Comparing 


0 

'Cl 
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The  overall  transfer  function  for  the  closed  loop  system  is 

■  -1  r 

XI  -  A,  0  Bs 

0  X[-Af  Bf 

=  C,{XI  -  A,)-^B,+Cfi\I-  Af)-^Bf 
=  rj(A)  +  r/(A)  (3. IS) 

There  are  three  sets  of  eigenvalues  for  the  asymptotic  transfer  function,  two  .s/oir  and 
one  fast.  The  first  set  of  stow  eigenvalues  is 

2i  =  {|a/, -/, +TA-f’A-2|  =  0} 

This  set  of  poles  is  asymptotically  uncontrollable  because  of  the  zero  element  in  B,- 

The  second  set  of  slow  eigenvalues  is  both  controllable  and  observable.  This  set 
contains  the  transmission  zeros  of  the  plant  augmented  with  the  minor  loop  rate  feedback 
through  the  measurement  matri.x,  M. 

Z2  =  I [a/„_;  -  I„_i  —  T Ai  +  T.-1i2A’2  'Fij  =  o| 

The  set  of  fast  eigenvalues  are  also  controllable  and  observable.  Note  that  eigenvalues 
obtained  from  .44  =  B2K1F2  and  F2B2K\  (which  is  used  below)  are  identical. 

23  =  {|A/,-/,  + =  0} 

.Vll  three  sets  of  poles  must  lie  in  the  open  left  half  s  plane  or  open  unit  circle  in  the 
:r  plane.  The  output  consists  of  the  fast  modes  containing  the  poles  and  a  set  of  slow 
modes  containing  the  poles  Z2. 
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3.3  DESIGN  PARAMETERS  FOR  THE  DISCRETE  FIXED-GAIN  PI  CONTROLLER 

In  order  to  achieve  a  diagonal  transfer  function,  r(A)  for  an  irregular  plant,  the 
measurement  matrix,  M  in  Figure  3.1  must  be  chosen  carefully.  .\I  should  be  formed  with 
as  few  non-zero  elements  as  possible  so  that 


F2  =  C2-\-MAx2  (3.19) 

has  full  rank  m  and  C2E2^  is  diagonal. 

The  following  steps  outline  the  procedure  to  determine  the  form  of  tlie  measiircineiit 
matrix  (3,  4); 


1.  Form  the  matrix 


5*  = 


clAi\Ax2 

.T  idn 


[3,20) 


where  m  is  the  number  of  control  inputs,  cj  is  the  fth  row  of  C\  and 


d,  =  mfn  :  cf .4|j,4i2  7^  0, j  =  0,  1,  ..,n-lj  (3.21) 

Equation  3.21  specifies  that  di  is  the  smallest  value  of  j  for  which  cJa\.^Ai2  #  0.  The 
permissible  values  of  j  are  0, 1, . .  .,n-  1.  If  all  the  values  of  j  result  in  cf.4ji.4i2  =  0, 
then  use  d,  =  n  -  1,  where  n  =  dimension  of  .4n. 

2.  Form  F2  =  C2  +  M A\2  using  a  general  form  of  the  matrix  M  =  The  elcmcnt.s 

m.j  appearing  in  F2  are  permitted  nonzero  values  only  if  D’  has  a  non-zero  element 
in  a  corresponding  position.  .411  other  elements  of  M  are  set  equal  to  zero. 

3.  Form  C2F2''  and  set  off-diagonal  terms  to  zero  if  possible  to  ensure  a  diagonal  matrix. 

4.  Ensure  that  F2  has  rank  m. 

For  the  CRC.4,  Equation  3.20  is  formed  using  the  parameters  from  .4CM  Entry. 
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0 

0 

1 

.0349 

0 

0 

0 

0 

0 

0 

0 

0 

0 

1 

0 

0 

1 

0 

0 

C2  - 

0 

0 

0 

0 

0 

0 

1 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

1 

0 

0 

0 

0 

1 

0 

0 

0 

111  = 

0 

0 

0 

II 

0 

0 

0 

1 

0.035 

-32.1804 

0 

-.0119 

-.0186 

-31.2 

15  0 

0 

0 

It  can  be  verified  that  the  form  of  M  does  not  change  for  the  different  fliglit  condi¬ 
tions  by  noting  the  same  location  of  non-zero  elements  in  .dn.  .412,  Cj.  and  C2  terms  of 
■dppendi.K  .4. 


Equation  3.20  yields  the  values  di  =  d2  =  f/3  =  d.^  =  d^  =  0.  Then  B  j  formed. 


-.0186  -31.235  0  0  0 

0  0  0  0  0 


0  10  0  0 
0  0  0  1  0.0349 

0  0  0  0  0 


(3.22) 


The  matrix  C2  in  Equation  3.19  has  a  dimension  5  x  5.  Since  the  dimension  of  .4i2  is  3  x 
5,  the  necessary  dimension  of  M  is  5  x  3.  Using  M  of  the  form 


m,i 

mi2 

"U3 

m2i 

17122 

777  23 

77332 

17141 

777.12 

^43 

mm 

77752 

^^53 

3-10 


the  matrix  Fi  of  Equation  3.19  is 


0349  -  .018677713 

7nii  - 

-  31.2357ni3 

0 

mu 

.03577712 

—  .0 1867n23 

777  21  ■ 

-  31.23577723 

1 

tJlni 

.035/7722 

—  .018677733 

7/731  - 

-  31.23577733 

0 

77732 

.0357773, 

-  .018677733 

77731  - 

-  31.23577733 

0 

777  3  2 

.035777.12 

-.018677753 

'»51  - 

-  31.23577753 

0 

mu 

1  4-  .035777 

The  assignable  elements  of  F2  are  permitted  non-zero  values  only  if  B’  has  a  nonzero 
element  in  the  same  position.  Thus,  it  is  required  that  m23  =  ^133  —  mi3  =  m33  =  m^]  = 
m.ii  =  rnsi  =  rui  i  =  ^22  =  "232  =  "^52  =  0.  This  leaves  an  /-h  matrix  of  the  form 


Fi  = 


.0349  -  .0186mi3  mu  —  31.235mi3  0  0  0 

0  0  10  0 

0  m3i  0  0  0 


0  Tn^2  .03577332 
0  0  1 


'3.2.5’ 


The  matrix  must  bo  diagonal,  therefore  the  third  element  in  the  first  I'ov.  >4' 

Equation  3.26,  rrzn  =  3l.235mi3-  It  i*  tilso  apparent  from  the  first  and  third  terms  of  the 
first  row  that  setting  77113  =  9  ensures  that  tire  matrix  product  is  diagonal  and  insensiti\'e 
to  variation  of  the  measurement  matrix  values.  With  non-zero  values  of  77731  ^itd  77?i2, 
rank  of  F2  is  m.  Since  the  measurement  matrix  implemented  in  the  computer  simulations 
has  only  two  non-zero  elements  77731  a-nd  77732,  a  reduced  5x2  M  matrix  is  used  since  the 
third  column  of  M  contains  only  zeros. 


0119 


0 

0 

0 


0 


.03^49  — 3 

0  i 


0 

0 

0 


0.3l9(mi|-.31.235mi..) 

-( .1)3-19— .0186  m  1 3)77131 

0 

0 

0 

0 


0  u 
0  0 
0  0 
0  0 
0  1 


(.’5.26) 
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The  proportional  path  gain  matrix,  in  Figure  3.1  is  selected  to  make  the  last 
transfer  function.  F/lA)  of  Equation  3.18  diagonal.  Therefore, 

=  E  (i.dTj 

where  E  is  a  diagonal  matrix  of  elements  a,.  So 

A'l  =  (Afi2)"^E 

and 

A  2  =  —  Aq  A  1  ( 3.21) ) 

Ao  is  the  ratio  of  integral  to  proportional  control. 

This  form  of  M  and  the  appropriate  choice  of  A'l  only  ensure  that  t!ie  fast  transfer 
function.  F/f  A )  of  Equation  3. 18  is  diagonal.  All  of  the  terms  of  Fsf  A )  must  also  be  diagonal 
to  ensure  an  overall  diagonal  transfer  function  FfA;. 

The  elements  of  the  E  matrix  are  selected  by  a  trial  and  error  process.  When  a 
diagonal  FfA)  has  been  achieved,  each  element  of  the  E  matrix,  a,  fine  tunes  the  ith 
input's  affect  on  the  ith  output. 

.3.4  DISCRETE  STEP  RESPONSE  PI  CONTROLLER  THEORY  AND  DESIGN 

Porter  has  developed  a  technique  that  can  be  used  to  design  a  control  system  for 
a  plant  when  the  mathematical  equations  are  not  known  (10.  8).  The  method  has  been 
developed  for  a  regular  and  asymptotically  stable  plants  (wdiere  all  open  loop  eigenvalues 
are  m  the  open  unit  circle  of  the  z  plane)  where  the  steady  state  transfer  function.  G(01, 
has  rank  m  ( no  transmission  zeros  at  the  origin)  and  A*,  the  plant  decoupling  matrix,  has 
full  rank  ni. 

cf.d'A  D 

A*  =  :  (3.30) 

m 
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rn  is  the  number  of  control  inputs,  cf  is  the  zth  row  of  C  and 


d,  =  min  ;  cj  .V  B  ^  O,^  =  0.  1 . n  -  l]  (:h:i! ) 

An  effort  has  been  made  in  Hammond's  thesis  to  extend  this  design  to  the  cu'^'  of 
an  irregular  and  unstable  plant.  The  extension  of  the  step  response  method  is  contiiuied 
in  this  thesis  effort. 

The  discretized  plant  Equations  3.4,  3.7  are  repeated  and  Equation  3.3  is  used  to 

form 

x(^-+l)T  =  <^i(kT)  + '^  u{fcT) 

:{k+l)T  =  :{kT)  +  Te(kT) 

=  :{kT)  +  Tr(kT)-Tw{kr) 

=  x(kT)  +  Tr(kT)  -TFx{kT) 

The  control  law  is 

u{kT)  =  r[A’ie(^'r)  +  h'2:{kT)]  (3.3'_bi 

.N’ote  that  Figure  3.1  is  applicable.if  the  gain  factor  is  changed  to  T  for  this  design. 
The  closed  loop  equation  can  now  be  expressed  as 

x{k  +  l)T  =  (<J>-Ttf'A'iF)x(fc7')  +  T'I'A'2^(/tT)  +  rTA'ir(/tr) 

:{k+i)T  =  -TFxikT)  +  z{kT)  +  Tr(kT) 

The  block  diagonalization  depicted  in  Equations  3.16  and  3,17  is  performed  and  tlio 
following  two  sets  of  eigenvalues  result.  The  sloiu  set  is  determined  from 

2,  =  i\\F  ~  U  +  TA\  =  0} 

and  the  fast  set  is 

22  =  {|aa  -At  t\a-^bK2\  =  o} 
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The  step  response  method  is  designed  around  two  plant  responses,  H(T)  and 
G(0)  =  H{oc).  The  response  with  a  step  input  at  the  lirst  sampling  time,  T,  is 

H{T)  =  -  [n}B 

Referencing  Equation  3.11,  H(T)  can  be  approximated  as  T  —  0  by 

H(T)  =  CA-K\TD  =  TCD 
Gis)  =  C(sIn-  A)-^B 

and  as  s  ^  0. 

G{0)=-CA-^B  (3.34) 

The  proportional  and  integral  path  gain  matricies  are  given  as  ( iO,  9.  S) 

h'l  =  (3.35) 

and 

A'2  =  G(0)-‘rn  (3.36) 

where  E  is  a  diagonal  matrix  as  defined  in  Equation  3.27  and  m  is  another  diagonal  matrix 
of  arbitrary  values  which  are  determined  by  trial  and  error  through  computer  simulations. 

Using  the  ,4,  B,  and  C  matricies  for  the  representative  flight  condition  of  Table  2.5, 
the  resulting  F’  and  G(0)  matricies  have  full  rank  of  m  =  5. 

If  the  plant  equations  are  known  in  advance,  the  approximations  of  Equations  3.33 
and  3.34  can  be  used.  Otherwise,  offline,  open  loop  tests  could  be  performed  on  the  plant 
to  determine  H(T)  and  G(0)  for  a  known  input. 

The  aircraft  equations 

X  =  ,4i  +  Bu 


and 


y  ~Cx 
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can  be  expressed  in  terms  of  ARM  A  coefRcients.  The  algorithm  developed  by  Velez  ill) 
was  "sed  to  generate  the  ARMA  (oefftcients  for  the  model. 

y{kr)  +  A\ii{k  —  l)T  ^  A  y(k  —  y  )T  =  B\u{k  —  IIT  +  •■•  +  B\  u.{k  —  )T  ( .3  ..37 ) 
It  can  be  seen  from  the  equation  above  that 

H(T) =  B, 

and  since,  y(0)  =  0. 

G(0)  =  (/+  .4i  +  ■•■  +  +  ■■■+  B^■) 

3.5  ADAPTIVE  STEP  RESPO.XSE  PI  CONTROLLER 

There  are  many  parameter  adaptive  algorithms  from  which  the  designer  can  clioixse. 
Hammond  chose  to  use  a  least  squares  method  since  “It  is  conceptually  simple  and  exhibits 
statistical  properties  that  are  as  good  as  those  of  maximum  likelihood  method  for  most 
practical  situations.”  (4:3-12) 

The  parameter  that  is  to  be  estimated  is  the  ARMA  coefficient  representation  of  the 
plant,  or  the  vector  9.  The  .A.RMA  model  representation  of  Equation  3.37  is  used  and  the 
.ARMA  output  Equations  2.6.,  2.7.  and  2.8  are  repeated  here. 

y{kT)  =  9'^  (kT)<t>(kT) 

where  for  the  reduced  second  order  ARMA  model  of  the  CRCA, 

9^{kT)  =  [aJ  bJ  bJ]^ 

oikT)'^  =  [->/{k  -  \  }T  -  /(A:  -  2)7  iJ[k  -  1)7  J{k  -  2)Tf 
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Hammond  used  the  following  equations  to  implement  the  least  squares  estimate  of  B 


e[k+  1) 
P(k+  1) 
't(k+  1) 


9{k)+ P{'-]T{k^i)^{k+i)[y(i-+i)-z^{hi  :)0{k]]  (;].:!8) 

P(k)  -  P{k)x(k  +  l)^(k  +  l)x^(k  +  l)Pik)  (3.30) 

_ 1 _ 

[Qf  +  x'^(k  +  l)P(k)x{k  +1)]  '  '*** 


where  i(A:  +  1)  is  Hammond’s  notation  in  the  adaptive  algorithm  for  d{k  +1).  P(k)  is  the 
covariance  matrix,  and  at  is  a  weighting  factor  =  1  for  this  effort. 

Calculation  of  P(k+1)  in  Equation  3.39  can  be  changed  by  the  introduction  of  a  for¬ 
getting  factor,  A  (Hammond  used  7  as  the  forgetting  factor,  but  that  notation  is  confusing 
since  there  is  another  term  ‘/{k  -f  1)  used  in  the  Equations  3.38  and  3.39). 


P{k+l)  =  \[Pik)  -  P{k}x{k  +  ly/ik  +  l)x^{k+ l)P{k)]  (3.41) 

A 

and  7(fc  +  1)  of  Equation  3.40  now  becomes 

7(^'+l)  [Aa( -I- x^(A: -f  l)P’(^')i'(/;  4- 1)] 

When  the  forgetting  factor  is  incorporated  in  simulations.  Equations  3.38,  3.41,  and 
3.42  are  used  to  update  the  changing  ARMA  coefficients. 

This  chapter  details  the  mathematical  development  of  Porter’s  discrete  PI  control 
laws;  discrete  PI,  step  response,  and  parameter  adaptive  based  upon  an  ARMA  model. 
Procedures  are  outlined  for  determination  of  the  controller  gain  matrices  A'l  and  A'2  based 
upon  known  plant  matrices  or  step  responses.  The  chapter  concludes  with  a  presentation  of 
the  recursive  least  squares  parameter  estimation  algorithm  to  be  used.  Chapter  4  cont.vins 
the  results  from  the  discrete  PI  controller. 
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[V.  DISCRETE  PI  CONTROLLER 


4 . 1  TRA  NSFER  F  UNCTION  A  -V.4  L  YSIS 

In  addition  to  the  plant  meeting  the  conditions  listed  in  Section  3.2,  the  steady  state 
control  surface  deflections  necessary  to  perform  each  maneuver  in  every  flight  condition 
(including  failures)  must  be  exaininea.  If  the  amount  of  control  surface  deflection  required 
in  steady  state  is  larger  than  the  position  limits  listed  in  Table  2.6,  then  the  maneuver 
cannot  be  performed  for  that  particular  flight  condition.  The  control  surface  rate  and 
position  limits  cannot  be  exceeded  at  any  time  in  the  simulations  or  the  responses  will 
become  nonlinear  when  limit  circuits  are  encountered.  The  steady  state  surface  deflections 
determined  for  each  combination  of  maneuvers  will  be  used  as  a  guide  to  determine  the 
feasibility  of  maneuver  performance. 

The  output  relationship  shown  below  is  used  to  determine  the  steady  state  surface 
deflections  (inputs)  required  to  perform  each  maneuver. 

y{s)  =  G(5)u(s) 

-Applying  the  final  value  theorem,  where  u„  is  the  magnitude  of  the  step  input  and  y„  is 
the  magnitude  of  the  desired  system  output, 

u„  =  G(0)-V„  (4.1) 

where  G(0)  is  defined  in  Equation  3.34,  which  is  repeated  below,  and  y,,  =  r(t)  for  a 
tracking  system. 

G(0)  =  -CA-^B 

The  five  maneuvers  and  their  associated  input  commands  are 

1.  pitch  rate  tracking  -  q=2°/sec  for  3  sec  and  zero  thereafter.  The  pitch  rate  command 

is  generated  by  ramping  6  to  6°  over  3  seconds. 


2.  coordinated  turn  -  45°  roll  angle  (d>)  and  yaw  rate  (r)  determined  from  Equation  4.2 

^  (S?) = 

where, 

g  =  the  gravitational  constant  {32.174 ftlsec^) 

V  =  the  forward  velocity  of  the  aircraft 
<i>  =  the  desired  bank  angle 

3.  sideslip  tracking  -  2°  of  sideslip  angle  (/3) 

4.  flat  turn  -  \°  j sec  yaw  rate  (r)  with  zero  roll  angle  {4>) 

5.  banked  turn  -  only  roll  angle  {<i>)  commanded 

Table  4.1  gives  the  corresponding  yaw  rate  (r)  for  each  of  the  flight  conditions  con¬ 
sidered  in  this  thesis  when  the  bank  angle  {4>)  is  equal  to  45  degrees. 

Table  4.1.  Vcmd  -  45°  Bank  Angle 


Flight  Condition 

Velocity 

ftfsec 

Gravity 
ftj  sec^ 

Yaw  Rate 
deg! sec 

.ACM  Entry 

895.0 

.ACM  Exit 

263.0 

iSB 

TF/TA 

1004.9 

32.17 

UEM 

The  output  vectors  for  each  maneuver,  y„,  consist  of  the  following  five  outputs  and 
the  units  are  expressed  in  ft/sec  for  the  velocity,  radians  for  the  angles,  and  rad/sec  for 
the  yaw  rate. 


V  =  forward  velocity 
0  =  sideslip  angle 
6  =  pitch  angle 
<f>  =  bank  angle 
r  =  yaw  rate 
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and  the  input  vector,  ti„,  consists  of  the  five  control  surface  inputs,  where  the  units  are 
expressed  in  degrees. 


8ci  =  left  canard 
=  right  canard 
Stei  =  left  trailing  edge 
Stgr  =  right  trailing  edge 
Srud  =  rudder 


The  steady  state  output  vector  for  a  45°  coordinated  turn  in  ACM  Entry  is: 


0.0000 

ft /sec 

0.0000 

rad 

0.1047 

rad 

0.0000 

rad 

0.0000 

rad/sec 

and  from  Equation  4.1,  the  steady  state  surface  deflections  are: 


-1.4411 

deg 

1.4411 

deg 

9.4767 

deg 

-0.4767 

deg 

-2.6545 

deg 

(4.3) 


(4.4) 


The  necessary  surface  deflections  to  perform  each  maneuver  in  all  flight  conditions 
considered  are  listed  in  Tables  4.2,  4.3,  4.4,  4.5,  4.6,  and  4.7.  The  G(0)  for  25%  loss  of 
rudder  in  ACM  Entry  is  not  invertible. 
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Table  4.2.  Steady  State  Control  Inputs  -  Pitch  Rate  Tracking 


Flight  Condition 

{deg) 

Scr 

{deg) 

^td 

{deg) 

^ter 

{deg) 

^rud 

ACMENTRY 

1.2100 

1.2100 

1.3836 

1.3836 

0.0000 

ACM30TL 

2.6967 

2.7687 

5.4502 

3.1944 

0.0000 

ACM50CL 

-16.3074 

-3.6573 

-4.6384 

-6.0154 

-3.6573 

ACM25RL 

N/A 

N/A 

N/A 

N/A 

N/A 

ACMEXIT 

-16.4536 

-14.4260 

-31.9885 

-28.6533 

-3.0725 

TFTA 

1.0384 

1.1469 

1.6532 

1.5927 

-0.0945 

TFTA30TL 

0.8025 

0.7998 

2.0189 

1.1727 

-0.0013 

TFTA50CL 

2.2060 

1.1090 

1.8015 

1.7633 

0.0060 

TFTA25RL 

-2.1867 

3.2095 

1.0801 

1.3309 

5.3962 

where, 

ACMENTRY  =  ACM  Entry  (nominal  flight  condition) 

ACM30TL  =  30%  loss  of  effectiveness  of  the  left  trailing  edge  -  ACM  Entry 
ACM50CL  =  50%  loss  of  effectiveness  of  the  left  canard  -  ACM  Entry 
ACM25RL  =  25%  loss  of  effectiveness  of  the  rudder  -  ACM  Entry 
ACMEXIT  =  ACM  Exit  flight  condition 
TFTA  =  TF/TA  flight  condition 

TFTA30TL  =  30%  loss  of  effectiveness  of  the  left  trailing  edge  -  TF/TA 
TFTA50CL  =  50%  loss  of  effectiveness  of  the  left  canard  -  TF/TA 
TFTA25RL  =  25%  loss  of  effectiveness  of  the  rudder  -  TF/TA 
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Table  4.3.  Steady  State  Control  Inputs  -  45°  Coordinated  Turn 


Flight  Condition 

6ci 

(deg) 

6cr 

(deg) 

hel 

(deg) 

^ter 

(deg) 

^rud 

id€g) 

.\CM  ENTRY 

-1.4411 

1.4411 

-2.6545 

ACM30TL 

-1.6473 

1.5466 

IB 

-2.8283 

ACM50CL 

0.8474 

1.6176 

HR  Bl 

wsm 

-1.2235 

ACM25RL 

N/A 

N/A 

N/A 

N/A 

N/A 

ACMEXIT 

46.8304 

42.2510 

100.3224 

90.4552 

7.7426 

TFTA 

-0.3664 

0.3664 

0.2181 

-0.2181 

-0.7977 

TFTA30TL 

-0.4190 

0.3684 

0.3435 

-0.2782 

-0.8498 

TFTA50CL 

-0.9488 

0.1031 

-0.1306 

-0.4729 

-0.6676 

TFTA25RL 

-1.9385 

1.9385 

0.2696 

-0.2696 

1.3886 

Table  4.4.  Steady  State  Control  Inputs  -  Sideslip  Tracking 


Flight  Condition 

^ci 

{deg) 

Scr 

(deg) 

(deg) 

^ter 

(deg) 

^rud 

(deg) 

ACMENTRY 

1.4671 

-1.4671 

0.2431 

BBSI 

4.2623 

ACM30TL 

0.8361 

-0.9397 

0.8124 

-3.4561 

ACM50CL 

-0.3633 

-0.6934 

0.1201 

-0.9094 

2.1692 

ACM25RL 

N/A 

N/A 

N/A 

N/A 

N/A 

ACMEXIT 

0.9444 

-0.8711 

1.6158 

-1.4159 

1.2350 

TFTA 

1.0973 

-1.0973 

0.2444 

-0.2444 

3.5678 

TFTA30TL 

0.9310 

-0.9946 

0.4317 

-0.3497 

3.3312 

TFTA50CL 

2.9287 

1.3314 

0.5312 

3.0699 

TFTA25RL 

3.5708 

0.2099 

-0.2099 

0.8878 
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Table  4.5.  Steady  State  Control  Inputs  -  Flat  Turn 


Flight  Condition 

(deg) 

Scr 

(deg) 

■Eiw 

^ter 

(deg) 

6rud 

(d^g) 

ACMENTRY 

9.3284 

-9.3284 

-3.0266 

3.0266 

15.7374 

ACM30TL 

11.0335 

-10.3688 

-5.2088 

3.8672 

17.4679 

ACM50CL 

-5.6600 

-10.8046 

-7.3813 

-4.9166 

6.6545 

ACM25RL 

N/A 

N/A 

N/A 

N/A 

N/A 

ACMEXIT 

7.4212 

-6.8452 

-5.1607 

6.7319 

15.2521 

TFTA 

2.9565 

-2.9565 

-1.6359 

1.6359 

5.0267 

TFTA30TL 

3.4114 

-3.0262 

-2.6129 

2.1166 

5.5084 

TFTA50CL 

7.7362 

-0.8407 

1.1951 

3.7252 

4.0173 

TFTA25RL 

13.2925 

-13.2925 

-1.8054 

1.8054 

-9.1315 

Table  4.6.  Steady  State  Control  Inputs  -  45°  Banked  Turn 
Flight  Condition  ^ter  ^rud 


{deg)  (deg)  (deg)  (deg)  (deg) 


.ACMENTRY 

-15.0300 

15.0300 

4.8857 

-4.8857 

-25.5798 

.4CM30TL 

-17.7202 

16.6512 

8.3768 

-6.2193 

-28.2744 

ACM50CL 

9.0926 

17.3571 

11.8667 

7.8892 

-10.9173 

ACM25RL 

N/A 

N/A 

N/A 

N/A 

N/A 

ACMEXIT 

36.0196 

52.2227 

107.8402 

80.6486 

-14.4756 

TFTA 

-4.2022 

4.2022 

2.3406 

-2.3406 

-7.3194 

TFTA30TL 

-4.8451 

4.2947 

3.7335 

-3.0243 

-7.9965 

TFTA50CL 

-10.9859 

1.1938 

-1.6811 

-5.3060 

-5.8797 

TFTA25RL 

-19.1844 

19.1844 

2.6120 

-2.6120 

13.2360 

Table  4.7.  Steady  State  Control  Inputs  -  15“  Banked  Turn 


Flight  Condition 

^cl 

(deg) 

6cr 

(deg) 

^tel 

(deg) 

(deg) 

^rud 

{deg) 

ACM ENTRY 

-10.0181 

10.0181 

-3.2565 

-17.0499 

ACM30TL 

-11.8112 

11.0987 

5.5835 

-4.1454 

-18.8460 

.4CM50CL 

6.0605 

11.5692 

7.9096 

5.2585 

-7.2738 

ACM25RL 

N/A 

N/A 

N/A 

N/A 

ACMEXIT 

24.0085 

34.8085 

71.8797 

53.7554 

-9.6486 

TFTA 

-2.8009 

2.8009 

1.5601 

-1.5601 

-4.8787 

TFTA30TL 

-3.2294 

2.8626 

2.4885 

-2.0158 

-5.3300 

TFTA50CL 

-7.3225 

0.7957 

-1.1205 

-3.5367 

-3.9190 

TFTA25RL 

-12.7871 

12.7871 

1.7410 

-1.7410 

8.8223 

Comparison  between  the  necessary  steady  state  control  surface  deflections  to  perform 
the  maneuvers  listed  in  Tables  4.2,  through  4.7  with  the  control  surface  limitations  listed  in 
Table  4.8  determine  which  maneuvers  can  possibly  be  performed  in  each  flight  condition. 


Table  4.8.  Control  Surface  Position  and  Rate  Limits 


Surface 

Position  Limit 
(deg) 

Rate  Limit 
(deg/sec) 

Canards 

-h60  -30 

±  100 

Trailing  Edges 

-1-30  -30 

±  100 

Rudder 

-f-20  -20 

±  100 

Comparison  between  Tables  4.2  -  4.7  and  Table  4.8  yields  the  fact  that  only  the 
sideslip  tracking  and  flat  turn  maneuvers  can  be  performed  in  ACM  Exit.  The  magnitude 
of  the  banked  turn  must  be  reduced  to  15“  in  ACM  Entry,  and  the  results  for  the  25% 
rudder  loss  conditions  are  questionable. 

4.3  TIME  RESPONSES 

The  time  responses  obtained  include  the  two  maneuvers  investigated  by  Hammond, 
pitch  rate  tracking  and  the  coordinated  turn,  plus  an  additional  three  maneuvers  listed 
in  Section  4.1.  .411  maneuvers  are  investigated  in  all  failure  conditions  possible,  including 
the  new  case  of  25%  rudder  loss.  A  single  set  of  design  parameters,  <T,,m3i,  77142,  Aq 
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rrom  Equations  3.28,  3.23,  and  3.29  are  used  for  all  tested  flight  conditions.  An  attempt  is 
made  to  obtain  a  universal  set  of  gains  that  will  stabilize  the  aircraft  in  all  combinations 
of  maneuvers  and  failure  situations.  The  values  of  the  design  parameters  are  chosen  to 
optimize  the  time  responses  and  to  try  to  keep  all  individual  controller  gain  values  below 
20.  The  gain  values  are  minimized  in  order  to  reduce  the  magnification  of  noise  signals. 


The  command  inputs  used  are  ramped  step  inputs  since  it  is  physically  unrealizable 
for  a  pilot  to  generate  a  pure  step  input.  The  following  approximation  has  been  used 
(4:5-14) 


ramp  time  = 


position  limit 


60° 


rate  limit  100°/sec 
where  the  largest  allowable  surface  deflection  for  the  canard  is  used.  A  ramp  time  of  0.5  sec 
is  used  for  all  inputs  to  simplify  the  MATRIX^  simulations.  In  addition  to  the  ramping, 
Hammond  also  inserted  a  prefilter  to  remove  high  frequency  spikes.  The  prefilter  is  selected 
to  limit  the  bandwidth  to  under  10  rad/sec  which  does  not  adversely  affecting  the  speed 
of  system  response. 


Prefilter  = 


10 


s  +  10 


.A.11  of  the  time  responses  collected  for  the  discrete  PI  design  used  the  following  values 
of  cr,,  m3i,  m42,  and  Aq  in  Equations  3.28  and  3.29: 


Oi 

=  0.2 

^31 

=  0.1 

<72 

=  0.05 

m42 

=  0.1 

<73 

=  0.5 

Aq 

=  0.3 

<74 

=  0.3 

<75 

=  0.5 

.Appendix  C  contains  the  K\  and  Ki  matrices  for  each  flight  condition  and  Table  4.9 
shows  the  stability  of  results  for  all  flight  conditions  using  the  design  parameters  listed 
above.  The  time  responses  have  been  generated  using  the  input  commands  shown  in 
Figures  4.1  through  4.5.  The  magnitudes  of  the  input  commands  (displayed  in  degrees) 
are  shown  for  ACM  Entry.  The  yaw  rate  differs  on  the  coordinated  turn  input  for  TFTA 
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(see  Table  4.1),  and  the  banked  turn  can  be  performed  at  45°  in  TFTA. 


Table  4.9.  Stability  Analysis  Using  discrete  PI  design  parameters 


Flight 

Condition 

Stable  With 
Universal  Gain 

Requires  Gain 
Change 

ACMENTRY 

X 

ACM30TL 

X  (I) 

ACM50CL 

X 

ACM25RL 

X 

ACMEXIT 

X 

TFTA 

X 

TFTA30TL 

X 

TFTA50CL 

X 

TFTA25RL 

X 

(1)  =  The  flat  turn  maneuver  is  not  stable  for  the  universal  gain. 


The  four  responses  in  Figures  4.1  through  4.5  are  defined  as  follows 

-  0  ...  9 

—  4>  ■-  r 


4-9 


theto  (deq) 


i-io 


L 

1 

1 

"  f 

“T - 

-  1 

-  1 

i 

J 

1 

1 

1 

7 

■ 

f  . 

_ _ _ 1 — i — 

»  '  • 

2  4  6  8  10  12  14 

seconds 

Figure  4..S.  Input  command  for  banked  turn  maneu 
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Thf  resulting  time  responses  for  simulations  in  ACM  Entry,  ACM  Exit,  and 
TFTA  follow  and  sample  surface  deflections  and  rates  are  included  only  for  the  first  ma¬ 
neuver.  pitch  rate  tracking.  Limiting  devices  are  included  in  the  simulation  to  keep  the 
surface  deflections  and  rates  within  the  limits  defined  in  Table  4.S.  Once  the  surface  de- 
flf?ction  and  rate  limits  are  e.xceeded,  the  limiting  devices  are  saturafed  and  tlie  siniuhit ioii 
becomes  nonlinear  and  this  normally  drives  the  time  responses  unstable. 

The  four  responses  on  each  of  the  ne.xt  plots  are  again  defined  as  follows 

—  3  ...  9 

- O  ■-  r 


Figure  4.6.  Di.crnte  Controller  pitch  rate  tracking  response  in  .‘\.CM  Entry 
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Figure  4.11.  Discrete  Controller  flat  turn  response  in  .AC.M  liiiti  v 


seconds 

Figure  4.12.  Discrete  Controller  banked  turn  response  in  .ACM  l.niiy 
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belc.thelo.phi.r  ( 


Figure  4.16.  Discrete  Controller  coordmated  turn  response  iii  I  FI  A 
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belo.thelo.phi.r  (deg.deq  per  sec 


Time  responses  for  a  selected  maneuver,  sideslip  tracking,  are  shown  for  all  of  rhe 
failure  conditions  in  ACM  Entry  and  TFTA.  The  simulations  were  performed  using 
tlie  lixed  A'l  and  AT  values  calculated  for  the  non-failed  ACM  Entry  and  TFTA  flight 
conditions  (see  Appendix  C).  From  Figures  4.22  and  4.25,  it  can  be  seen  that  the  side.-lip 
tracking  task  cannot  be  accomplished  for  the  25%  rudder  loss  in  ACAI  Entry  or  TFT.A. 


Figure  4.20.  Discrete  Controller  sideslip  maneuver  response  in  .\C.M30TL 
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Figure  4.21.  Discrete  Controller  sideslip  maneuver  response  in  AC.MiOCL 


Figure  4.22.  Discrete  Controller  sideslip  maneuver  response  in  ACM2.5RL 


beta, Iheto, phi, r  (deg, deg  per  sec 


seconds 


Figure  4.25.  Discrete  Controller  sideslip  maneuver  response  in  TFT.-\25R.L 
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4.3  FREQ UE.\C Y  .4 .V.4 L  YSIS 

Conventional  frequency  analysis  of  a  MIMO  system  is  difficult,  but  is  made  easier 
when  fhe  outputs  are  relatively  decoupled,  as  is  the  case  with  this  design  method.  Since  the 
outputs  are  decoupled,  only  the  diagonal  transfer  fuction  pairs  are  considered,  i.e, 

•For  s  atically  unstable  aircraft,  the  military  specification  MIL  -F  -9490D,  specifies  the 
gain  margin  to  be  ±  6.0  dB  and  the  phase  margin  to  be  ±  45  degrees.  For  these  types  of 
aircraft  the  system  may  have  a  low  and  high  frequency  gain  margin.  .4t  the  low  frequency 
crossover  frequency  of  -180  degrees  the  gain  margin  should  be  smaller  than  -6.0  dB.  .A.t  the 
high  frequency  crossover  of  -180  degrees  the  gain  margin  should  be  greater  than  4-6.0  dB 
f  1: 4-17).”  The  Nichols  plots  of  0/6cmd  in  Figures  4.28  and  4.32,  illustrate  the  high  and  low 
gain  margin  characteristics  of  the  unstable  CRCA. 

The  air-'raft/controller  model  simulated  in  MATRIX^  is  a  hybrid  model  consisting  of 
a  co'.itinuous  plant  with  a  discrete  controller.  In  order  to  conduct  a  frequency  analysis,  an 
equi%-alent  discrete  model  of  the  overall  system  must  be  generated.  This  is  accomplished  by 
treating  the  controller  as  a  continuous  unit  and  then  transforming  the  “continuous  system” 
to  the  w/  plane  via  a  Hofmann  transformation.  The  Hofmann  transformation  used  in  this 
thesis  was  obtained  from  Hammond’s  computer  files  (4:4-29). 

Flw')  = 

Gill')  = 

H{w')  = 

D(w')  = 


■Vccording  to  Houpis  and  Lamont  (5),  quantities  in  the  wf  plane  can  be  handled  like 
continuous  quantities  and  the  z  plane  roots  can  be  found  from  the  following  equation. 


Selected  Nichols  and  Bode  diagrams  are  included  for  both  the  healthy  ACM  Entry 
and  TFTA  flight  conditions.  Figures  4.26  through  4.33  show  the  open  loop  frequency 
characteristics  of  the  discrete  PI  controller  design.  The  closed  loop  Bode  plots  show  the 
bandwidth  of  the  system  with  the  control  law  incorporated  for  the  two  selected  outputs, 
in  Figures  4.34  through  4.37.  The  bandwidth  characteristics  of  the  healthy  aircraft  design 
are  included  in  Table  4.10  and  the  tvf  plane  roots  are  listed  in  Table  4.11. 

lible  4.10.  Bandwidth  of  discrete  PI  design 


output 

ACM  Entry 

TFTA 

i; 

7.22 

0.30 

9 

9.77 

4> 

10.23 

7.74 

r 

9.54 

8.90 

Table  4.11.  wi  plane  roots  for  .\CM  Entry  and  TFTA 


ACM  ENTRY 

TFTA 

-0.1119  ±j  0.1536 

•0.2266  ±j  0.1380 

-0.2960 

-0.3017 

-0.3014 

-0.3020 

-0.3055 

-0.3074 

-0,4947 

-0.7228 

-2.0425 

-3.6890 

-8.2710 

-9.9482  (5) 

-9.9482  (5) 

-11.0842  ±j  5.1379 

-7.4002  ±]  10.4555 

-13.1893 

-14.7533  ±j  13.2824 

13.1976 

-38.8237 

-29.2678  ±;  35.2494 

-40.7383  ±j  25.8551 

-16.4291  ±;  .54.7333 

-36.6469  ±j  33.5002 

-57.4723 

-50.2105  ±j  29.3661 

-51.5347  ±j  30.4394 

-52.7656  ±]  30.7213 

-61.4256 

-.58.1497  ±]  31..3455 

-57.9341  ±j  31.3275 
-14.8236  ±]  74.5446 
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Figure  4.26.  Discrete  FM  .Mchols  plot  of  in  .AC.M  l  .iiirv 


!■  igiire  4.27.  Discrete  PI  open  loop  Pode  plot  of  J/.j. 
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Figure  4.28.  Discrete  PI  Nichols  plot  of  in  AC.M  liniry 
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f  igure  1.29.  Discrete  PI  i»pen  loop  Bode  [)lot  of  m  .\(  M  f.nM  v 
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Figure  -1.30.  Discrete  PI  Nichols  plot  of  !■>  1  1  TA 
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Figure  4.31.  Discrete  PI  open  loop  Bode  plot  of  in  I  I  l  A 
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Figure  4.34.  Discrete  PI  closed  loop  Bode  plot  of  i^l 3cmd  in  ACM  Entry 
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Figure  4.3.5.  Discrete  PI  closed  loop  Bode  plot  of  d/Ocmd  in  .ACM  I'.nti  v 
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Figure  4.36.  Discrete  PI  closed  loop  Bode  plot  of  i3/,3cmd  in  TFT.4 
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Figure  4.37.  Discrete  PI  closed  loop  Bode  plot  of  9/9cmd  in  1  F  I.V 
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The  phase  and  gain  margins  for  ACM  Entry  and  TFTA  are  listed  in  Table  4.12. 


Table  4.12.  Gain  and  Phase  Margins  -  Discrete  PI  Design 


ACM  Entry 


Transfer 

Function 

Gain  \ 
(dl 

largin 

B) 

UJc 

(rad /sec) 

Phase 

Margin 

(deg) 

U4, 

(rad /sec) 

Low 

High 

Low 

High 

V 

(2) 

14.26 

(2) 

70.62 

74.64 

8.07 

P 

(2) 

48.54 

(2) 

50.54 

(1) 

(1) 

9 

-25.22 

12.65 

2.48 

64.44 

39.66 

21.51 

<t> 

(2) 

17.07 

(2) 

64.34 

44.80 

14.47 

r 

(2) 

8.60 
_ 1 

(2) 

50.84 

52.45 

20.10 

TF/TA 


Transfer 

Function 

Gain  Margin 
(dB) 

(rad /sec) 

Phase 

Margin 

(deg) 

(rad /sec) 

Low 

High 

Low 

High 

V 

-39.18 

16.51 

0.24 

50.40 

70.61 

7.21 

P 

(2) 

48.54 

(2) 

50.69 

(1) 

(1) 

9 

-65.15 

6.63 

0.05 

43.22 

43.75 

22.38 

4> 

(2) 

18.66 

(2) 

70.68 

57.91 

13.80 

r 

(2) 

8.73 

(2) 

51.56 

126.45 

0.55 

180.03 

0.74 

54.45 

20.35 

(1)  =  Response  is  always  less  than  0  dB 

(2)  =  No  low  frequency  gain  margin  (pha.se  >  -180  degrees) 

This  chapter  contains  a  listing  of  the  steady  state  control  surface  deflections  to  per¬ 
form  each  maneuver  in  the  chosen  flight  conditions.  The  design  parameters  and  time 
responses  for  the  discrete  PI  controller  are  presented  for  non  failure  and  single  failure 
conditions.  Frequency  analysis  of  the  discrete  design  is  presented  for  the  non  failed  cases 
ACM  Entry  and  TFTA.  Chapter  5  contains  the  time  responses  for  the  step  response 
controller  and  adaptive  algorithm. 
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STEP  RESPONSE  METHOD  AND  ADAPTIVE  RESULTS 


5. 1  S TEP  RESPONSE  METHOD 

The  time  responses  collected  for  the  step  response  method  use  the  same  inputs  dis¬ 
played  in  Figures  4.1  through  4.5  and  the  following  values  of  cr,,  77131,  ^7142,  and  777,  in 
Equations  3.35  and  3.36. 


=  0.1 

77ri 

0.2 

02 

=  0.015 

77^2 

= 

0.2 

=  0.01 

7tr3 

= 

0.2 

O4 

=  0.05 

7^4 

= 

0.2 

=  0.075 

= 

0.2 

mai 

=  0.1 

77142 

= 

0.1 

The  w'  plane  roots  are  calculated  for  the  step  response  method  in  the  same  fashion 
as  discussed  for  the  discrete  PI  controller  in  Section  4.3  and  are  included  in  Appendix  C 
along  with  the  fixed  gain  values  of  K\  and  K2  used  in  the  simulations  for  the  no  failure 
cases  ACM  Entry  and  TFTA.  The  computer  simulations  are  conducted  with  a  gain  of 
1.  The  gain  factor  in  the  control  law  of  Equation  3.32,  T,  is  already  incorporated  in  the 
gain  matrices  listed  in  Appendix  C. 

Table  5.1  shows  the  stability  results  for  all  the  possible  maneuvers  (see  Section  4.1) 
in  the  various  flight  conditions  for  the  step  response  method  design  parameters.  The  same 
time  responses  are  shown  here  as  have  been  presented  for  the  discrete  PI  controller.  The 
four  responses  on  each  plot  are  distinguished  according  to  the  same  legend  used  in  Chapter 
4. 

-  0  ...  9 

—  ,p  —  r 
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Table  5.1.  Stability  Anal 


ysis  Using  step  response  design  parameters 


Flight 

Condition 

Stable  With 
Universal  Gain 

Requires  Gain 
Change 

ACMENTRY 

X 

ACM30TL 

X 

ACM50CL 

X 

ACM25RL 

X 

ACMEXIT 

X  (1) 

TFTA 

X 

TFTA30TL 

X 

TFTA50CL 

X  (2) 

TFTA25RL 

X 

(1)  =  The  responses  are  slightly  divergent  for  ACM  Exit  over  the  20  second  simulation 
time. 

(2)  =  All  results  have  a  distinct  periodic  oscillation  and  slowly  diverge  over  the  simulation 
duration. 
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Figure  5.7.  Step  response  m«  thod  flat  turn  response  in  .AC.M  K.xii 
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Figure  -5  S  Step  response  method  pitch  rate  tracking  response  in  Tl' IW 
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Figure  5,15.  Step  response  method  sideslip  tracking  response  in  .\CM25RL 


Figure  5.16.  Step  response  method  sideslip  tracking  response  m  FI-  F.ARD'l  1. 
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5.2  ADAPTIVE  CONTROLLER 


When  a  single  set  of  design  parameters  for  a  fixed-gain  design  do  not  adequately 
stabilize  and  control  the  aircraft,  either  gain  scheduling  or  an  adaptive  algorithm  can  be 
used.  Hammond’s  adaptive  algorithm,  which  is  built  on  the  step  response  method,  is 
examined. 


Since  the  aircraft  plant  can  be  represented  by  a  discrete  difference  equation 
for  a  given  sampling  time  T,  control  law  gains  are  calculated  by  manipulating 
the  coefficients  of  this  difference  equation.  By  measuring  input  and  output 
data,  on-line,  the  coefficients  that  constitute  the  difference  equation  itself  can 
be  updated  continuously  and  new  PI  controller  gains  calculated  each  time  the 
measurements  are  taken.  Predicting  the  new  plant  parameters  allows  gain 
calculations  to  be  accurate  and  appropriate  for  the  new  flight  condition  or 
configuration.  (4:6-1) 


.A.  recursive  least  squares  (RLS)  algorithm  is  used  to  update  the  ARMA  coefficients  each 
sampling  period  which  are  in  turn  used  to  determine  the  controller  gain  matrices. 

The  RLS  algorithm  is  written  in  FORTRAN  executable  code  and  accessed  each 
sample  period  in  the  MATRIXjf  simulation.  The  algorithm  is  initialized  with  a  covariance 
matrix  P{0)  of  I  and  the  matrix  containing  the  past  output  and  input  values,  i(0)=  0. 
.ARMA  coefficients  for  either  the  non-failed  ACM  Entry  or  TFTA  flight  conditions  are 
used  for  the  first  two  sampling  periods,  0.05  seconds,  to  ensure  that  the  x{k  +  1)  matrix 
is  fully  populated  before  the  ARMA  coefficients  in  0{k  -|-  1)  are  estimated.  The  rate  of 
calculation  of  the  varying  ARMA  coefficients  is  controlled  by  the  choice  of  forgetting  factor, 
A.  The  applicable  RLS  Equations  3.38,  3.41,  and  3.42  are  repeated  he'e  with  the  weighting 

eik)  +  P{k)xik  +  l)7(it-|-  l)[y{k  +  1)  -  x^{k  +  l)eik)] 

i[F(A:)  -  Pik)xik  +  l)7(fc  +  l)i^(A:  +  l)P(fc)] 

_ 1 _ 

[Aof  -h  -I-  l)P{k)x{k  -I-  1)] 

adaptive  algorithm  is  listed  in  Appendix  D. 


factor,  Q(  =  l. 

5(it+l)  = 
P{k  +  l)  = 
l(k+l)  = 

The  source  code  for  the 
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T!ie  adaptive  controller  is  initially  evaluated  by  analyzing  the  results  of  all  of  the 
liiaiieiivers  performed  in  ACM  Entry  and  TFTA  with  no  failures,  hhe  loi'gc'l  1 1 tig  lactc)! . 
A  is  set  to  one  and  the  ARMA  coefficiLiits  are  monitored  at  the  onset  and  upon  completion 
of  the  maneuvers.  The  .ARM.A  coefficients  do  not  change  when  the  adaptive  .ilgorithiii 
i.'  luii.  which  is  the  e.xpected  result.  The  outputs  are  e.xpected  to  closely  lesemhle  tho.-^f 
oi)i ained  from  the  stop  response  method,  since  the  adaptive  simulations  start  with  the  same 
.\RMA  plant  models  used  in  the  step  response  simulations.  The  following  time  losponse-,. 
however,  do  not  show  ‘he  desired  decoupling  of  outputs  that  has  been  obtained  from  hot  !i 
the  discrete  PI  and  step  response  design  methods. 

The  four  responses  on  each  of  the  ne.xt  plots  are  delinod  as  follows 
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Figure  .5.19.  Adaptive  pitch  rate  tracking  response  in  .\('.\I  Fntry 
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Figure  5.20.  .adaptive  coordinated  turn  response  in  ACM  F.iiiiy 
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Figure  5.21.  .Adaptive  sideslip  tracking  response  in  .ACM  Fniry 
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Figure  5.22.  .Adaptive  flat  turn  response  in  .ACM  Entry 
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Figure  5. 2.1.  .Adaptive  banked  turn  response  in  .ACM  Entry 
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Figure  5.26.  Adaptive  sideslip  tracking  response  in  TF  IW 
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Figure  5.27.  .Adaptive  flat  turn  response  in  FF  12V 
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Figure  5.28.  .Adaptive  banked  turn  respon.se  in  TF  I'.A 
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Adaptive  simulations  are  now  performed  with  a  single  failure  occuiiing  ai  ciMi't 
thf  lioginning  of  the  simulation  or  at  two  seconds.  The  valiu's  tliat  are  lead 

into  M AT RIXx  to  initialize  the  0  vector  are  for  non-failed  cases.  The  simulations  are 
performed  for  a  variety  of  A  values  in  the  single  failure  cases  of  dO'T  left  trailing  edge 
loss  and  509c  left  canard  loss  for  the  sideslip  tracking  and  flat  turn  maneiners.  Only  iln' 
responses  generated  in  the  failure  condition  .4CM30TL  are  presented  since  i  he  oi  h'T  ( ,!-e.^ 
■piinluccd  'dvergent  responses.  For  the  two  maneuvers  tested  in  the  iwi,  ddfi'i'eni  f.idiin' 
oii-.'is,  the  responses  obtained  from  a  A  =.90  are  the  most  stabh'.  I  he  re-pon''  ~ 
i  nn  Ic.ir  30  seconds  for  the  case  of  failure  introduction  at  two  seconds  in  oiiiej-  to  ailow  liic 
responses  to  settle. 

Ihe  four  responses  on  each  plots  are  defined  according  to  the  fo’lo-.s'ing  h'uend. 
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at  1  =  0 


Deto.tneto.pni.r  (aeg.oeq  per  sec) 


iiiure  5.30.  .Adaptive  flat  turn  response  in  ACM.30TEL  wiili  A  = 


beto, theta, phi, r  (deg, deg  per  sec) 


4 
3 
2 
1 
0 

-1 

0  3  6  9  12  15  18  21  24  27  30 

seconds 

Fi>>,ure  5.32.  Adaptive  flat  turn  response  in  ACM30TEL  with  A  =  .90  and  failure  at  1=2 
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The  adaptive  controller  is  initially  e\-aluatcd  by  analyzing  the  r<-suUs  of  all  of  the 
:::aiieuvers  performed  in  ACM  Entry  and  TFTA  with  no  failures.  The  loigeiiiug  facKn. 
\  -el  to  one  and  the  ARMA  coeiricienls  are  monitored  at  the  onset  and  ii[ion  complei ioii 
Ilf  the  maneuvers.  The  ARMA  coefficients  do  not  change  wlion  the  adapiivi'  algoinhni 
i-  iiin,  which  is  the  expected  result.  The  outputs  are  expected  to  closely  re.-emidi'  tho-' 
o!)tained  from  the  step  response  method,  since  the  adaptive  simulations  start  with  tiie  s.-imi- 
.Ali.MA  plant  models  used  in  the  step  response  simulations.  The  following  time  response-, 
h'iwever.  do  not  show  the  desired  decoupling  of  outputs  that  has  been  obtained  fi  om  boi  h 
fin'  discrete  PI  and  step  response  design  methods. 

The  lour  responses  on  each  of  the  next  plots  are  deliiied  a-  follov.s 
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Figure  .5.19.  .Adaptive  pitch  rate  tracking  response  in  .AC.M  Entry 
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Ail.iptivc  simulations  are  now  performed  with  a  single  lailtire  nccuriiue  .ii  i 
!i"  iioginning  of  the  simulation  or  at  two  secon.ds.  The  ARMA  valuo  Miat  aic 
!i!o  M  A  T R[ X X  to  initialize  the  B  vector  are  for  non-failed  cases.  The  siiuulati(ui 
ii'il'ornied  for  a  variety  of  A  values  in  the  single  failure  cases  of  'M)X.  h'ft  trailing 
M^s  and  505t  left  canard  loss  for  the  sideslip  tracking  and  flat  turn  iiianeii\fr.s.  On! 

. opcrn.ses  generated  in  the  failure  condition  ACM30TL  are  [rreseiiteil  -inci'  A'le  m  la  r 
irodiicod  divergent  responses.  For  the  two  maneuvers  tested  in  tlu'  t'.'.o  d.:!'  :!  ::'  I. 
iii''  is,  the  responses  obtained  from  a  A  =.90  are  the  most  stable.  1  he  re-pda-e- 
!ii,  for  30  seconds  for  the  case  of  failure  introduction  at  two  seconds  in  order  to  .i!hv 
■* s- ponses  to  settle. 

The  four  responses  (tti  each  plots  are  defined  according  to  the  followine  l.'o.aid. 
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f  igure  .‘j.gO.  .\d;iptive  sideslip  tracking  response  in  .'\('M.30  1  f.l.  with  \-  'lO  a'ld  I 
at  t=i) 
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I  ii;:::*'  '>.'.V2.  A(ia[iti\'e  flat  turn  r<'S[)nnso  in  ACM-'fO'rKI  with  A  =  .9()  tmd  failure  at  t=2 
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This  chapter  contains  the  design  parameters  and  time  responses  for  the  step  response 
PI  controller  and  results  of  the  adaptive  algorithm.  The  results  of  all  work  in  this  thesis 
are  analyzed  in  Chapter  6. 
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VI.  CONCLUSIONS  AND  RECOMMEND.ATIONS 
6.1  SUMMARY 

6.1.1  DISCRETE  PI  CONTROLLER  The  time  domain  and  frequency  results  ob¬ 
tained  for  the  discrete  PI  controller,  presented  m  Chapter  4,  show  that  the  outputs  are 
indeed  highly  decoupled.  The  velocity  results  are  not  plotted  in  the  same  figures  as  the 
other  four  outputs  because  the  180°/7r  radian  conversion  fa^'^or  is  not  applicable.  The 
velocity  responses  are  monitored  and  never  peaked  above  2  ft/sec,  which  is  well  within 
an  acceptable  range  of  variation.  The  time  responses  for  the  non  failed  flight  conditions 
show  that  the  design  is  very  robust.  The  bandwidth  of  the  system  was  increased  by  a 
factor  of  100  over  the  previous  design  and  is  now  just  under  10  radfsec  for  most  of  the 
five  input/output  pairs.  The  phase  and  gain  margins  specified  by  MIL  -  F  -  9490D  are 
met  in  almost  all  inrtances.  The  goal  of  determining  one  set  of  design  parameters  that 
can  be  used  to  determine  fixed  controller  gains  for  all  three  flight  conditions  are  met  and 
robust  results  are  obtained  for  ACM  Entry,  TFTA,  and  ACM  Exit.  An  effort  is  made 
to  limit  the  values  of  the  individual  elements  in  the  Ki  and  K2  matrices.  Examination 
of  .Appendix  C  demonstrates  that  the  ACM  Entry  and  TFTA  gains  are  in  a  reasonable 
range  and  only  a  few  of  the  elements  in  the  ACM  Exit  gain  matrices  are  larger  than 
desired.  The  extension  of  the  fixed-gain  design  to  the  Ccises  of  single  failure  conditions 
is  fairly  successful  in  that  for  three  of  the  six  fcdlure  conditions,  the  aircraft  does  remain 
stable  and  perform  the  maneuver  with  some  degree  of  flying  qualities  intact.  However, 
some  type  of  gadn  scheduling  scheme  or  adaptive  control  is  deemed  necessary  by  the  fart 
that  three  failure  conditions  yield  unstable  responses. 

6.1.2  STEP  RESPONSE  METHOD  PI  The  step  response  method  PI  controller 
time  domain  results  show  the  satme  degree  of  stability  and  decoupling  between  outputs, 
using  one  set  of  design  parameters,  as  are  obtained  from  the  discrete  PI  controller.  The 
step  response  method  shows  two  improvements  over  the  discrete  PI  controller  design;  the 
3  undershoot  in  the  flat  turn  is  reduced  by  approximately  one  third  in  ACM  Entry 
and  TFTA,  and  the  flat  turn  can  now  be  successfully  performed  in  ACM30TL.  The 
goal  of  obtaining  individual  gain  elements  of  less  than  20  is  not  possible  in  this  method 
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and  the  step  response  controller  can  only  stabilize  two  of  the  six  single  failure  conditions. 
Examination  of  the  u;’  roots  listed  in  Appendix  C  aJso  shows  that  all  roots  are  not  in  the 
LHP. 


6.1.3  .ADAPTIVE  CONTROLLER  The  results  of  the  adaptive  controller  are  not 
all  a.s  expected.  The  simulations  performed  for  the  no  failure  cases  show  much  larger 
0  outputs  in  the  sideslip  tracking  responses  and  slightly  larger  4>  responses  in  the  flat 
turn  maneuver  than  are  in  the  step  response  method  time  responses.  This  result  occurs 
even  though  the  ARMA  coefficients  are  invariant.  This  characteristic  of  the  4)  output 
undoubtedly  carries  over  into  the  single  control  surface  failure  simulations.  The  steady 
state  values  of  the  adaptive  outputs  are,  however,  very  close  to  the  steady  state  step 
response  outputs. 

Hammond  had  worked  with  the  coordinated  turn  in  ACM30TL,  so  two  different 
maneuvers  in  two  single  failure  conditions  are  attempted.  However,  the  responses  can  not 
be  stabilized  in  ACM50CL,  TFTA30TL,  and  TFTA50CL.  It  is  interesting  to  note  that 
it  is  the  9  response  which  is  consistently  divergent  in  the  TFTA  failure  cases.  The  steady 
state  outputs  for  the  two  sets  of  fculure  onsets  are  reasonably  decoupled  at  the  completion 
of  the  simulation. 

6. 2  RECOMMENDA  TIONS  FOR  FUTURE  RESEA  RCH 

There  are  many  areas  that  can  be  examined  using  Porter’s  control  law  techniques  on 
the  CRCA.  The  first  of  which  are  areas  of  immediate  interest  and  deal  with  the  parameter 
adaptive  algorithm. 

1.  Determine  why  there  are  differences  between  the  sideslip  tracking  and  flat  turn  ma¬ 
neuvers  obtained  from  the  step  response  PI  method  and  adaptive  simulations. 

2.  Stabilize  the  adaptive  simulations  in  the  two  failure  conditions  attempted  in  this 
thesis. 

.3.  Assess  the  effects  noise  in  the  output  measurements  have  on  the  fixed-gain  and  adap¬ 
tive  designs. 
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4.  Use  non-linear  aircraft  equations  to  perform  simulations  with  the  design  parameters 
determined  from  the  linear  point  designs  and  compare  the  results. 

5.  Investigate  the  proportional  plus  integral  plus  derivative  (PID)  control  law  on  the 
same  set  of  maneuvers  and  failure  conditions  to  remove  the  slow  modes  of  the  time 
responses. 

6.  Compare  the  results  obtained  from  the  PI  design  techniques  to  those  obtained  using 
the  PID  control  law. 
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Appendix  A.  AIRCRAFT  STATE  SPACE  MODELS 


The  aircraft  models  used  in  the  MATRIX^  simulation  are  listed  on  tlie  followimr 
pages. 
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Table  A.l.  ACM  Entry  Matrices  -  No  Failures 


.0000 

.0000 

.0000 

.0000 

1.0000 

.0000 

.0000 

.1)000  ■ 

,0000 

.0000 

.0000 

.0000 

.0000 

.0000 

1.0000 

u.!.:o 

-32.1804 

.0000 

-.0119 

-.0186 

-31.2350 

.0000 

.0000 

0000 

-1.0634 

,0000 

-.0324 

-1.0634 

894  4548 

.0000 

.0000 

0000 

.0000 

.0000 

.0000 

.0069 

-.6015 

.0000 

.0000 

0000 

.0000 

.0360 

.0000 

.0000 

.0000 

-.0929 

.0349  - 

OC  \ 

.0000 

.0000 

.0000 

.0000 

.0000  - 

27.8066 

-2,0376 

.4913 

.0000 

.0000 

.0000 

.0000 

.0000 

2.4582 

-.0241  - 

.4377  . 

.0000 

0000 

OOOU 

.0000 

.0000  .0000 

,0000 

.0000 

.0000 

1 

0000 

0000 

0000 

0000 

.0000  .0000 

.0000 

.0000 

.0000 

.0411 

0411 

.1322 

0866 

.1322  .0866 

.1018 

.1018 

,0000 

-.3163  - 

3163  - 

.9597  - 

.6194  - 

.9597  -.6194 

-1.0183 

-1.0183 

.0000 

.1014 

1014  - 

.0284  - 

.0215  - 

.0284  -.0215 

-.0200 

-.0200 

.0000 

.0003  - 

0003  - 

.0002  - 

.0001 

.0002  ,0001 

-.0001 

.0001 

,0006 

.0762  - 

0762 

.2219 

2011  - 

.2219  -.2011 

.1109 

-.1109 

,1144 

.0486  - 

0486 

.0029 

0021  - 

.0029  -.0021 

.0021 

-.0021 

-.0544 

■  .0000 
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1.0000 

.0349 
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.0000 
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.0000 
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B  Matrix  for  Five  Control  Surfaces 


.0000 

.0000 

.0000 

.0000 

.0000 

.0000 

.0000 

.0000 

.0411 

.0411 

.3206 

.3206 

-.3163 

-.3163 

-2.5974 

-2.5974 

.1014 

.1014 

-.0699 

-.0699 

.0003 

-.0003 

-.0004 

.0004 

.0762 

-.0762 

.5339 

-.5339 

.0486 

-.0486 

.0071 

-.0071 
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Table  A. 2.  ACM  Entry  Matrices  -  30%  Loss  Of  Effectiveness  Left  Trailing  Edge 


,A 


B 


.0000 

.0000 

.0000 

.0000 

1.0000 

.0000 

.0000 

.0000  ■ 

.0000 

.0000 

.0000 

.0000 

.0000 

.0000 

1.0000 

.0450 

-.32.1420 

.0000 

-.0050 

.0550 

-39.9760 

.0000 

.0000 

.0000 

-1.4370 

.0000 

-.0240 

-1.0280 

894.1070 

.0000 

.0000 

.0000 

.0000 

.0000 

.0000 

.0070 

-.6920 

.0000 

.0000 

.0000 

.0000 

.0360 

.0000 

.0000 

.0000 

-.0990 

.0450 

-.9990 

.0000 

.0000 

.0000 

.0000 

.0000 

-31.8340 

-2.1380 
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.0000 

.0000 

.0000 

.0000 

.0000 

2.6670 

-.0310 

-.4420  . 
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.0000 

.0000 

.0000 

.0000 

.0000 

.0000 
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.0000 

.0000 

.0000 

.0000 

.0000 

.0000 

.0000 

.0000 

.0000 

.0520 

.0520 
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.0780 

.1040 

.0670 
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.0800 

.0000 

-.3300 

-.3300 

.0000 

-.6220 

-.9650 

-.6220 

-.7350 

-.7350 

.0000 
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.1020 

.0000 

-.0210 

-.0280 

-.0210 

-.0200 

-.0200 

.0000 

.0000 

.0000 

.0000 

.0000 

.0000 

.0000 

.0000 

.0000 

.0010 

.0800 

-.0760 

.0000 

.2010 

—.2220 

-.2010 
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-.1110 

.1150 

.0490 

-.0480 

.0000 
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-.0040 

-.0030 
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-.0030 
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■  .0000 

.0000 

1.0000 

.0410 

.0000 

.0000 

.0000 

.0000 

.0000 

.0000 

.0000 

.0000 

.0000 

1.0000 

.0000 

.0000 
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1.0000 

.0000 

.0000 

.0000 

.0000 

.0000 

.0000 

.0000 

.0000 
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.0000 
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.0000 
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B  Matrix  for  Five  Control  Surfaces 
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.0000 

.0000 

.0000 

.0000 

.0000 

.0000 

.0000 

.0000 

.0000 

.0000 

.0520 

.0520 

.1710 

.2510 

.0000 

-.3303 

-.3300 

-1.3570 

-2.3220 

.0000 

.1020 

.1020 

-.0410 

-.0690 

.0000 

.0000 

.0000 

.0000 

.0000 

.0010 

.0800 

-.0760 

.3120 

—  .5340 

.1150 

.0490 

-.0480 
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-.0100 

-.0550 
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Table  A. 3.  ACM  Entrv  Matrices  -  50%  Loss  Of  Effectiveness  Left  Canard 


.0000 

.0000 

.0000 

.0000 

1.0000 

.0000 

.0000 

.0000  ■ 

.0000 

.0000 

.0000 

.0000 

.0000 

.0000 

1.0000 

.0390 

-.32.1500 

.0000 

-.0080 

.0580 

—36.6690 

.0000 

.0000 

.0000 

-1.2460 

.0000 

-.0290 

-1.0450 

894.3280 

.0000 

.0000 

.0000 

.0000 

.0000 

.0000 

.0040 

-.6460 

.0000 

.0000 

.0000 

.0000 

.0360 

.0000 

.0000 

.0000 

-.0820 

.0390 

-.9990 

.0000 

.0000 

.0000 

.0000 

.0000 

-23.8380 

-2.0700 

.4990 

.0000 

.0000 

.0000 

.0000 

.0000 

2.4580 

-.0270 

-.4520  . 

.0000 

.0000 

.0000 

.0000  .0000 

0000 

0000 

.0000 

.0000  • 

.0000 

.0000 

.0000 

.0000  .0000 

0000 

0000 

.0000 

.0000 

-.1560 

-.1150 

.1320 

.0840  .1320 

0840 

1000 
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.0000 

-.1150 

-.2250 

-.9610 

-.6210  -.9610  - 

6210  - 
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-.7320 

.0000 

.0360 

.0720 

-.0230 

-.0220  -.0280  - 

0220  - 

0200 
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.0000 

.0000 

-.0010 

.0000 

.0000  .0000 
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0000 
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.0350 

-.0690 
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.2010  -.2220  - 
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.1140 

.0240 
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-.0020 

—  .0550  . 

■  .0000 
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1.0000  .0349  .0000 

.0000 

.0000 

.0000  ■ 

.0000 

.0000 

.0000  .0000  .0000 

1.0000 

.0000 

.0000 

= 

1.0000 

.0000 

.0000  .0000  .0000 

.0000 

.0000 

,0000 

.0000 

1.0000 

.0000  .0000  .0000 

.0000 

.0000 

.0000 

.0000 

.0000 

.0000  .0000  .0000 

.0000 
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B  Matrix  for  Five  Control  Surfaces 


.0000  ,0000 

.0000 

.0000 

.0000 

.0000  .0000 

.0000 

.0000 

.0000 

-.1560  -.1150 

.3160 

.3160 

.0000 

-.1150  -.2250 

-2.3140 

-2.3140 

.0000 

.0360  .0720 

-.0700 

-.0700 

.0000 

.0000  -.0010 

.0000 

.0000 

.0010 

.0350  -.0690 

.5340 

-.5340 

.1140 

.0240  -.0490 
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-.0080 

-.0550 
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Table  A. 4.  ACM  Entry  Matrices  -  25%  Loss  Of  Effectiveness  Rudder 


uooo 

.0000 

.uOuo 

.UOOO 

.0000 

.0000 

.0000 

.0000 

-32.1540 

.0000 

-.0080 

.0540 

-1.1340 

.0000 

-.0320 

-1.0580 

.0000 

.0000 

.0000 

.0070 

.0000 

.0360 

.0000 

.0000 

.0000 

.0000 

.0000 

.0000 

.0000 

.0000 

.0000 

.0000 

.0000 

.0000 

.0000 

.0000 

.0000 

.0000 

.0000 

.0000 

.0410 

.0410 

.1320 

.0870 

-3170 

-.3170 

-.9600 

-.6200 

1020 

.1020 

-.0280 

-.0220 

.0000 

.0000 

.0000 

.0000 

.0760 

-.0760 

2220 

.2010 

.0000 

.0000 

.0000 

.0000 

1.0000 

.0000 

.0000 

.0000  ■ 

.0000 

.0000 

1.0000 

.0350 

-31.5470 

.0000 

.0000 

.0000 

894.4400 

.0000 

.0000 

.0000 

-.6740 

.0000 

.0000 

.0000 

.0000 

-.0920 

.0350 

-.9990 

.0000 

-27.8550 

-2.0330 

.4600 

.0000 

2.0520 

-.0270 

-.4190 

0000 

.0000 

.0000 

.0000 

.0000  ■ 

0000 

.0000 

.0000 

.0000 

.0000 

1320 

.0870 

.1020 

.1020 

.0000 

9600 

-.6200 

-.7330 

-.7330 

.0000 

0280 

-.0220 

-.020U 

-.0200 

.0000 

0000 

.0000 

.0000 

.0000 

.0000 

2220 

-.2010 

.1110 

-.1110 

.0870 

0000 

.0000 

-2.0520 

.0270 

.1190  . 

.0000 

.0000 

1.0000 

.0349 

.0000 

.0000 

.0000 

,0000 

.0000 

.0000 

.0000 

.0000 

.0000 

1.0000 

.0000 

.0000 

1.0000 

.0000 

.0000 

.0000 

.0000 

.0000 

.0000 

.0000 

.0000 

1.0000 

.0000 

.0000 

.0000 

.0000 

.0000 

.0000 

.0000 

.0000 

.0000 

.0000 

.0000 

.0000 

.0000 

1.0000 

B  Matrix  for  Five  Control  Surfaces 


.0000 

.0000 

.0000 

.0000 

.0000  1 

.0000 

.0000 

.0000 

.0000 

.0000 

.0410 

.0410 

.3210 

.3210 

.0000 

-.3170  - 

-.3170 

-2.3130 

-2.3130 

.0000 

.1020 

.1020 

-.0700 

-.0700 

.0000 

.0000 

.0000 

.0000 

.0000 

.0000 

.0760  - 

-.0760 

.5340 

-.5340 

.0870 

.0000 

.0000 

-2.0520 

.0270 

.4190 
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Table  A. 5.  ACM  Exit  Matrices  -  No  Failures 


.0000 

.0000 

.0000 

.0000 

1.0000 

.0000 

.0000 

.0000 

.0000 

.0000 

.0000 

.0000 

.0000 

.0000 

1.0000 

.5184 

-28.5877 

.0000 

-.2762 

-.0283 

-136.2190 

.0000 

.0000 

.0000 

-4.9465 

26.9479 

.1330 

-.6753 

262.7934 

.0000 

.0000 

.0000 

.0000 

.0000 

-.0018 

.0035 

-.6511 

.0000 

.0000 

.0000 

.0472 

.0322 

.0000 

.0000 

.0000 

-.0245 

.4602 

-.8878 

.0000 

.0000 

.0000 

.0000 

.0000  - 

-7.7280 

-1.4530 

.9687 

.0000 

.0000 

.0000 

.0000 

.0000 

-.0889 

-.0543 

.0456 

.0000 

0000 

.0000 

.0000 

.0000 

.0000 

.0000 

.0000 

.0000 

.0000 

.0000 

.0000 

.0000 

.0000 

.0000 

.0000 

.0000 

.0000 

-.0253  - 

.0281 

.0206 

.0131 

.0206 

.0131 

.0154 

.0154 

.0000 

-.1512  - 

.1471  - 

.1302  - 

.0848  - 

.1302 

-.0848 

-.0997 

-.0997 

.0000 

.0161  - 

.0156  - 

.0024  - 

.0018  - 

.0024 

-.0018 

-.0016 

-.0016 

.0000 

.0007  - 

.0007  - 

.0001  - 

.0001 

.0001 

.0001 

-.0001 

.0001 

.0004 

.0113  - 

.0337 

.0270 

.0240  - 

.0270 

-.0240 

.0135 

-.0135 

.0277 

.0133  - 

.0121  - 

.0005  - 

.0004 

.0005 

.0004 

-.0004 

.0004 

-.0130 

.0000  .0000  1.0000  .4600  .0000  .0000  .0000  .0000  ■ 

.0000  .0000  .0000  .0000  .0000  1.0000  .0000  .0000 

1.0000  .0000  .0000  .0000  .0000  .0000  0000  .0000 

0000  1.0000  .0000  .0000  .0000  .0000  .0000  .UOOO 

.0000  .0000  .0000  .0000  .0000  .0000  .0000  1.0000 

B  Matrix  for  Five  Control  Surfaces 


.0000 

.0000 

.0000 

.0000 

.0000  ■ 

.0000 

.0000 

.0000 

.0000 

.0000 

-.0253 

-.0281 

.0491 

.0491 

.0000 

-.1512 

-.1471  - 

.3147 

-.3147 

.0000 

.0161 

.0156  - 

.0058 

-.0058 

.0000 

.0007 

-.0007  - 

.0003 

.0003 

.0004 

.0113 

-.0337 

.0645 

-.0645 

.0277 

.0133 

-.0121  - 

.0013 

.U013 

-.0130  _ 
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Table  A. 6.  TF/TA  Matrices  -  No  Failures 


.0000 

,0000 

,0000 

.0000 

1,0000 

.0000 

.0000 

.OOuO  ■ 

.0000 

.0000 

.0000 

.0000 

.0000 

,0000 

l.UOOO 

.0153 

-3'2.1961 

.0000 

—  .0355 

.0357 

-15.6105 

.0000 

.0000 

.0000 

-..3002 

,0000 

-.0071 

-3.2056 

1004.8788 

.0000 

.0000 

.0000 

.0000 

.0000 

-.0003 

.0202 

-1.6773 

.0000 

.0000 

.0000 

.0000 

.0320 

.0000 

.0000 

.0000 

-.2538 

.0155 

-  9999 

.0000 

.0000 

.0000 

.0000 

.0000 

-66.9300 

-5.4612 

1.0349 

.0000 

.0000 

.0000 

.0000 

.0000 

8.2821 

-.0299 

-1.2709 

,0000 

.0000 

.0000 

.0000 

.1000 

.0000 

.0000 

,0000 

QoOO 

.0000 

,0000 

.0000 

0000 

.0000 

.0000 

.0000 

,0000 

.0000 

-  3284  - 

-.3284 

,2548 

.1679 

.2548 

.1079 

.1975 

,197 

5 

onno 

-6788  - 

-.6788 

-1.7517 

■1.1313 

-1.7517  -1 

.1313 

-1.3351  - 

1,3351 

.0000 

.2387 

,238' 

-  0534 

-.0406 

-.0534 

0406 

-.0372 

-.037 

,0000 

0012  - 

-.0012 

-.0001 

-.0001 

.0001 

.0001 

-UOOI 

,0003 

.0018 

2336  - 

-  2336 

4200 

3737 

-.4200 

.3737 

.2100 

-.2100 

,3737 

1590  - 

-  1590 
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Appendix  B.  ARMA  MODEL  GESERATIOS 


III  IXTRODUCTION 

This  appendix  contains  a  suininary  of  the  ARMA  model  conversion  techniip:''  de. 
'.•doped  by  Bokor  and  Keviezky  and  the  MATRIX^-  macro  used  to  accomplish  li.is 
hirriiation  for  all  the  plant  matrices  (1,  11).  Detailed  rnathemaiical  dev i  lopm.en: 
pi'K-edure  is  contained  in  reference  (1). 

B.J  ARMA  MODEL 

I  h.c  rnetliod  used  to  generate  the  .-VRMA  model  representation  from  a  state-space 
moviel  representation  is  based  upon  using  constructibility  invariants.  This  melhud  elimi- 
natos  some  of  the  problems  associated  with  the  observability  matrix  tocliniipie  ■viimi  the 
^tafe  transition  matrix  is  singular.  Recall  that  the  .state-s[)ace  represeutaiioii  of  the  rim- 


t'nuous  time  system  is  expressed  by 

X 

=  Ax(t)  +  Bu{t) 

y<t) 

=  Cx{t) 

vvliere 

.1  =  the  continuous  time  plant  matrix  (n  x  n) 

B  =  the  continuous  time  control  matrix  (n  x  m)  with  the  rank  or  B  = 
('  -  the  continuous  time  output  matrix  (/x  n) 
r  -  the  state  variable  vector  with  estates 
■i  =  the  control  input  vector  with  m  inputs 
V  =  the  output  vector  with  /  outputs 


The  state-space  matrices  of  hquation  B.l  and  Equation  B.2  are  discretized  for  a 


given  sampling  time  T  and  expressed  in  the  discrete  time  domain  by 


i{T)  =  Fx(T)  +  GuiT} 

(B.;5 

y(T)  =  Hx{T) 

(B.-l 

where 

F  =  the  discrete  plant  matrix  (n  x  n) 

G  —  the  discrete  control  matrix  (n  x  m)  with  the  rank  =  m 
H  =  the  discrete  output  matrix  (/  x  n) 

X  =  the  state  variable  vector  with  n  states 
u  =  the  control  input  vector  with  m  inputs 
y  =  the  output  vector  with  /  outputs 

For  this  technique,  the  discrete  time  representation  given  in  Equations  B.3  and  B.  l 
must  have  no  poles  at  the  origin,  i.e.  the  inverse  matrix  exists.  The  constructibility 
matrix  Co  is  defined  as 


and  must  have  rank  Cq  =  n. 


Co 


II F-^ 
HF-^ 
HF-^ 


IIF-n 


( n.'j) 


new  matrix,  TT ,  is  formed  from  the  linearly  independent  rows  of  in  E(pia- 
tinn  B.5  (where  each  block  H  F~'  is  /  x  n).  The  first  row  of  the  new  matrix  i.^^  the  first 
row  of  Co-  The  second  row  of  T  is  the  next  linearly  independent  row  of  Co-  The  process 
continues  until  the  matrix  TT  has  full  rank  n. 
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The  transformation  matrix  T  is  formed  by  rearranging  the  rows  of  matrix  TT  into 
tlie  following  form 


hiF-^ 

T  =  : 

hiF-^ 

hiF-'^' 

where  hi  is  the  ith  row  of  the  H  matrix. 


For  I  <  i  <  I,  the  ith  constructibility  index  V,  is  defined  as  the  smallest  positive  integer 
such  that  h,F~^'~^  is  a  linear  combination  of  the  rows  before  it.  Then  the  constructibility 
indices  satisfy  the  relation 

Vi  +  V2  +  •  •  ■  +  Vi  =  n  (B.7) 


and  are  arranged  in  descending  order 


Vi>V2>V3>...Vi  (B.3) 

If  Equation  B.S  is  not  satisfied,  permutate  the  output  matrix,  H ,  to  satisfy  the  relationship 
of  the  constructibility  indices. 

Then  the  following  matrices  are  defined: 

R  =  T"*  (B.9) 

//  =  HR  (B.IO) 

F  =  R-^FR  (BID 

G  =  R-^G  (B.r2) 
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A  reduced  order  ARMA  model  is  used  where,  instead  of  having  n  coelficionts,  only 
.V  =  n/rn  (rounded  to  the  largest  integer)  are  needed.  The  .-1,  coefficients  of  the  .-MIM.A. 
model  are  calculated  from  the  relationship 


A,  =  -HSpXk-i}  (B.l.d) 

wliere  i  =  1 . V,  (A;  -  i)  is  a  delay  operator,  and 


1 

I 

0 

0 

{k-2) 

0 

0 

{k-l\) 

0 

0 

0 

1) 

0 

0 

{k-2) 

0 

0 

(k-V2) 

0 

0 

0 

...  {k-D 

0 

0 

...  {k-2) 

0 

1 

0 

...  _ 

.A.S  an  example,  to  find  the  .4i  coefficients  using  Equation  B.13  and  Equation  B.M, 
the  Sp^{k)  matrix  elements  associated  with  the  (k-1)  terms  in  Equation  B.14  would  be  set 
equal  to  1  and  aJl  other  matrix  elements  equal  to  zero.  The  A2  coefficients  would  require 
Sp^ik)  matrix  elements  associated  with  the  (k-2)  terms  in  Equation  B.14  be  set  equal  to  1 
and  all  other  matrix  elements  equal  to  zero,  etc. 

The  Bi  coefficients  of  the  .ARMA  model  are  calculated  from  the  relationship 

B,  =  HS\{k-i)G  (B.IO) 
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where. 


where 


N 


1=1 


5  =  block  diagonal  [5i,  Sj, . . . ,  5;] 


and  Sj  is  the  block-diagonal  Toeplitz  matrix  given  by 


(B.16) 


(B.17) 


for  j  =  1, . . . ,  / 


5,  = 


0 

1 


(B.18) 


B.3  ARM  A  IMPLEMENTATION 

The  implementation  of  the  theory  presented  in  the  previous  section  is  executed  in 
M.A.TRIXx  using  a  modification  of  a  macro  developed  by  Velez  in  his  thesis  (11).  For  this 
design, 


V;  =  2 

V2  =  2 

Vj  =  2 

K,  =  1 

Ks  =  1 


B-5 


^Pl 


1  0  0  0  0 
0  0  0  0  0 
0  10  0  0 
0  0  0  0  0 
0  0  10  0 
0  0  0  0  0 
0  0  0  1  0 
0  0  0  0  1 


'p-i 


0  0  0  0  0 
1  0  0  0  0 
0  0  0  0  0 
0  10  0  0 
0  0  0  0  0 
0  0  10  0 
0  0  0  0  0 
0  0  0  0  0 


1  0  0  0  0  0  0  0 
0  1  0  0  0  0  0  0 
0  0  1  0  0  0  0  0 
0  0  0  1  0  0  0  0 
0  0  0  0  1  0  0  0 
0  0  0  0  0  1  0  0 
0  0  0  0  0  0  1  0 
0  0  0  0  0  0  0  1 


= 


(B.19) 


(B.20) 


(B.21) 
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00000000 
1  0  0  0  0  0  0  0 

00000000 
0  0  1  0  0  0  0  0 

5,3=  =5*  (B.22 

00000000 

0  0  0  0  1  0  0  0 

00000000 

00000000 

B.3.1  \fATRIXx  AR.\fA  Macro  The  following  MATRlX.v  inacio  lisling  can  be  dl 
rectly  typed  into  MATRIX a'  and  implemented  by  using  the  command  ••EXEC( 'filename')" 


//ENSURE  THAT  A/C  MATRICES  FOR  DESIRED  FLIGHT  CONDITION  HAVE  BEEN 
//LOADED  IN  MATRIXX. 

//THIS  PROGRAM  WAS  MODIFIED  FROM  VELEZ’S  AND  HAMMOND’S  THESES  AND  IS 
//WRITTEN  FOR  AN  8  STATE  PLANT  WITH  5  OUTPUTS 

N=8;  //DEFINE  THE  NUMBER  OF  STATES 

L=5:  //INDICATE  NUMBER  OF  OUTPUTS 

//THE  C  MATRIX  FOR  THIS  MODEL  HAS 

CDUM=C(Cl  2345],:);  //lO  OUTS,  5  OF  WHICH  ARE  COMMANDED 

//FORM  THE  S  PLANE  SYSTEM  MATRIX 

SDUM=<A,B:CDUM,0*EYE(5)>:  //CONTINUOUS  SYSTEM  MATRIX 

SD=DISC(SDUM,N, .025) ;  //DISCRETIZES  THE  SYSTEM  MATRIX 

//TSAMP=.025 

[F  G  H  D]=SPLIT(SD,N) ;  //SEPARATES  DISCRETE  MATRICES 
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//FORM  CONSTRUCTABILITY  MATRIX 

J=N-L:  //J+1  IS  THE  LARGEST  POSSIBLE  VALUE 

//OF  INDEX  VI 

//THE  REMAINING  SECTIONS  OF  THE  CODE 
//ARE  WRITTEN  FOR  THE  UNDER  THE 
//ASSUMPTION  THAT  J>0. 

CO=<H*INV(F)>:  //THIS  LINE  FORMS  THE  COMPLETE 

//CONSTRUCTABILITY  MATRIX  IF  J=0. 


FOR  r=2: J+1,  .  .  . 

co=<co:H*iNV(F)**i>: . . . 

END; 


//ONLY  J+L+1  ROWS  OF  CO  ARE  NEEDED 
//TO  SELECT  BASIS  VECTORS.  THIS 
//FORMS  THE  REMAINING  BLOCKS  OF  CO 


RCO=RANK(CO) 


//MAKE  SURE  THAT  CO  HAS  RANK  N 


/  /  ****t^**t*-tt************^*********1f**********************!^**ir******** 

//FORM  TEMPORARY  TRANSFORMATION  MATRIX 


KEEP=1: 

CR=0: 

TT=<CO(l, :)>: 
TEMP=TT; 


//TRACKS  WHICH  ROWS  OF  CO  WHICH  GO  IN 
//TEMP  INITIALIZES  COUNTER 
//INITIALIZES  TEMPORARY  VECTOR  1ST  HOW 
//INITIALIZES  TEMPORARY  VECTOR 


//LOOP  FOR  GENERATING 
//TEMPORARY  TRANSFORMATION  MATRIX 
FOR  I=2:J*L+1,  .  .  .  //LOOK  THROUGH  THE  ROWS  OF  CO  FOR 

TEMP=<TT:CO (I,  .  //THE  1ST  N  LINEARLY  INDEPENDENT  ROWS. 

RN=RANK(TEMP):...  //THIS  ENSURES  THAT  EACH  ROW  OF  THE 

IF  RN=I-CR  THEN  TT=TEMP ; . . .//H  MATRIX  IS  REPRESENTED. 
KEEPERS»<KEEP;I>; . . .  //"KEEP"  LISTS  ROWS  OF  CO  WHICH  ARE 

KEEP=KEEPERS : . . .  //INCLUDED  IN  "TEMP" . 

ELSEIF  RN<I-CR  THEN  TEMP=TT;... 

CR=CR+1: . . . 

END ...  . 

END; 


RN 


//MAKE  SURE  THAT  RANK  OF  TEMP  IS  N 
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/  /  *****************************************************:^i^7^:^^i:^^m*^^^^i^^ 

//CALCULATE  INDICES 


V=0*ONES(L,1) : 


//INITIALIZES  INDEX  MATRIX 


FOR  I=1:L, . . . 

K=KEEP(I.l)  :  .  .  . 

IF  K=I  THEN  V(I,1)=1;...  //ROWS  1-L  OF  CO  MUST  BE  LINEARLY 
END; .  .  .  //INDEPENDENT  FOR  THE  L  INDICES  TO 

END;  //EACH  EQUAL  AT  LEAST  1. 


FOR  I=L+1:N,...  //THIS  LOOP  IS  NOT  NEEDED  FOR  J=0 

K=KEEP(I,1) : . . .  //CASE. 

IF  K<(J-1)*L+1  THEN  V(K-L. 1)=V(K-L, 1)+1 : • • • 

ELSEIF  K<J’fL+l  THEN  V(K-(J-1)*L)=V(K-(J-1)*L)  +  1 ;  .  .  . 

END;  .  .  . 


IF  K=J*L+1  THEN  V( 1 , 1) =V( 1 , 1)  +  1 ;  .  .  . 

END;...  //PLEASE  NOTE  THAT  ADDITIONAL  LOOPS 

END;  //MUST  BE  ADDED  FOR  CASE  WHERE  J>3 

//OR  SUBTRACTED  WHEN  J<3. 


KEEP  //ensure  THAT  ALL  L  INDICES  AT  LEAST 

//ARE  EQUAL  1,  THEIR  SUM  IS  N,  AND 

V  //THEY  ARRANGED  IN  DESCENDING  ORDER. 

//OTHERWISE,  THE  H  MATRIX  MUST  BE 
//RESTRUCTURED . 

//THE  FOLLOWING  PORTION  OF  THIS  MACRO  WORKS  FOR  THE  5  INDICES  HAVING 
//VALUES  OF  2, 2, 2, 1,1  RESPECTIVELY.  OTHERWISE,  THE  "T"  MATRIX  AND 
//SP’S  AND  SQ  MAY  HAVE  A  DIFFERENT  FORM.  REFER  TO  BOKOR  AND 
//KEVICZKY  ARTICLE  FOR  SOECIFIC  GUIDANCE. 


T=<TEMP(1, :) ;TEMP(6, :) ;TEMP(2, :) ;TEMP(7. :) ;TEMP(3, :) ;TEMP(8, :) ; . . . 
TEMP(4, ;) ;TEMP(5, :)> ;  //ROWS  OF  TEMP  ARE  ARRANGED  IN  ORDER 

//OF  CALCULATED  INDICES 


//T=<H(1 , : )*F**1 ; . . . ;H(1 , : )*F**V1; | . . . ; | ; 1 

//  H(L. :)*F**1; . . . ;H(L, :)*F**VL> 
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R=INV(T) ; 

HBAR=H*R: 

FBAR=INV(R)*F*R; 

GBAR=INV(R)*G; 

SQ1=EYE(8);  //SQ1=S**0 


SQ2=<0  0000000;. 
10000000;. 
00000000;. 
00100000;. 
00000000;. 
00001000;. 
00000000;. 
0000000  0>; 


//SQi  IS  AN  NXN  BLOCK  DIAGONAL  MATRIX 
//WITH  L  ENTRIES  DIMENSIONED  ViXVi 
//HAVING  I'S  ON  THE  LOWER  SUB- 
//DIAGONAL  (A  VERSION  OF  A  TOETLITZ 
//MATRIX).  SQi=S**(i-l) 

//S=BLOCK  DIAG0NAL{T0EPLITZ(V1XV1)  , 

//  . . . .TOEPLITZ(VLXVL)>. 


//THE  NUMBER  OF  ARMA  COEFFICIENTS  NEEDED  FOR  A  REDUCED  ORDER  MODEL 
//IS  APPROXIMATELY  K=N/L  (ROUND  UP  TO  NEXT  INTEGER  FOR  ANY  REMAINDER) 


//FOR  THE  N=8  L=5  CASE,  2  ARMA  COEFF’S  AND  2  SP  MATRICES  ARE  NEEDED. 


//FOR  THE  K=1  CASE 

SP1=<1  0  0  0  0; . . . 

0  0  0  0  0; . . . 
01000;... 
00000;... 
00100;... 
00000;... 
00010;... 

0  0  0  0  1>; 

//AND  FOR  THE  K=2  CASE 

SP2=<0  0  C  0  0; . . . 
10000;... 

0  0  0  0  0;  .  .  . 
01000;... 

0  0  0  0  0;  .  .  . 
00100;... 

0  0  0  0  0; . . . 

0  0  0  0  0>; 


//THE  NXL  SP  MATRICES  HAVE  THE 
//SAME  FORM.  THE  FIRST  COLUMN  IS 
//CK-1;...;K-V1;0*(V2,1); . . .;0*(VL,1)] 
//THE  SECOND  COL  ENTRIES  START  ON  THE 
//ROW  UNDER  THE  FIRST  COL  ENTRIES 
//C0*(V1,1) ;K-1, . . . ;K-V2;0*(V3,1) ; .  .  .] 
//THIS  TREND  IS  CONTINUED  FOR  ALL  L 
//COLUMNS. 


//THE  SP  MATRICES  HAVE  A  TIME 
//SHIFT  OPERATOR  (K-i). 


//NOW  ALL  THE  PRELIMINARY  VARIABLES  HAVE  BEEN  DETERMINED  AND 
//CALCULATION  OF  THE  ARMA  COEFFICIENTS  CAN  PROCEED. 
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B1ARMA=HBAR*SQ1*GBAR: 
B2ARMA=HBAR*SQ2*GBAR : 


//Bl  COEFFICIENT 
//B2  COEFFICIENT 


A1ARMA=-HBAR*SP1: 

A2ARMA=-HBAR*SP2: 


//Al  COEFFICIENT 
//A2  COEFFICIENT 


G0=INV(EYE(5)+A1ARMA+A2ARMA)*(B1ARMA+B2ARMA) ; 


CLEAR  CDUM  SDUM 
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TO  rec;r]:s:sive  moving  average  models 


Table  B.l.  ACM  Entry  ARMA  Model  -  No  Failures 


-I.9766D  -  00 
O.OOOOD  -  00 
3.0638D  -  03 
O.OOOOD  -  00 
O.OOOOD-00 

9,76660-01 
0.00000  -  00 
-3.06270  -  03 
0,00000  -  00 
O.OOOOO-OO 

7,34840  -  0-1 
-6,78970  -  06 
3.15340  —  05 
2.40020  -  05 
1.20780  -  03 

-7.71840  -  04 
-2.16950  -  05 
3.37060  -  05 
-6.09300  -  04 
1.39380  -  05 


0.00000  -  00 
-1.98770  -  00 
0.00000  -  00 
-27.82600  -  00 
5.75490  -01 

0.00000  -00 
9.89830  -  01 
0.00000  -00 
27.79940  -  00 
-6.35930  -  01 

7.3484  0  -  04 
6.78970  -  06 
3.15340  -  05 
-2-40020  -  05 
-1.20780  -  03 

-7.71840  -  04 
2.16950  -  05 
3.37060  -  05 
6,09300  -  04 
-1.39380  -  05 


2.02000  -  00 
0.00000  -  00 
-1.98570  -  00 
0.00000  -  00 
0.00000  -  00 

-2.00150  -00 
0.00000  -  00 
9.88140  -  01 
0.00000  -  00 
0.00000  -  00 

5  80750  —  u3 
-6.38970  -  06 
-2.17880  -  05 
1.64150  -  04 
l.'VloO  04 

-5.54160  -  03 
1.34840  -  05 
-4.13800  -  06 
3.78700  -  04 
-8.66300  -  06 


0.00000  -  00 
-7.07180  -  06 
0.00000  -  00 
-9.74960  -  01 
-6.00380  -  04 

0.00000  -  00 
0.00000  -  00 
0.00000  -  00 
0.00000  -  00 
O.OOOOO  -  00 

5.80750  —  o.i 
6.3897  0  -  06 
-2.17880  -  05 
-1.61150  -  04 
-1.72360  -  04 

-5.54160  -  03 
-1.3484  0  -  05 
-4.13800  -  06 
-3.7870  0  -  04 
8,66300  -  06 


O.OOOOO  -  00 
—  8.7077  O  —  05 
0.00000  -  00 
-7.00190  -  01 
-9.73260  -  01 

O.OOOOO  -  00 
O.OOOOO  -  00 
0.0000  0  -  00 
0.0001)0  -  00 
O.OOOOO  -  00 

O.UOOOO  -  00 
.!. 31300  -  0.5 
0,00000  -  00 
.1.44220  -  05 
-1.3526  0  -  03 

o.OOoOO  -  00 
.3.2708  0  -  06 
0.00000  -  00 
9.1860  0  -  05 
-2.10130  -  06 


=  [.4i(l, l),Ai(l, 2), ...,.4i(5,5),. 42(1,1) . 42(5,5), 

fli(l.  1) . i?i(5,5), /?2(1M), - Z?2(5.5)] 
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Table  B.2.  ACM  Entry  ARMA  Model  -  30%  Loss  of  Effectiveness  Left  Trailing  Edge 


-1.9278T)  -  00 

U.OOOOO  -  00 

7.20580 

-  01  0.00000  - 

-  00 

o.ouooO 

-  Ill) 

O.OOOOD  -  00 

-1.98790  -  00 

0.00000 

-  00  -1.44820  - 

-  05 

1.13160 

-  04 

-1.3075D  -  02 

0.00000  -  00 

-2.03340 

-  00  0.00000  - 

-  00 

0.00000 

-  00 

O.OOOOD  -  00 

-21.45700  -  00 

0.00000 

-  00  -9.80680  - 

-  01 

-5,40420 

-  01 

O.OOOOD-00 

5.62540  -  . 

0.00000 

-  00  -5.95870  - 

-  04 

-9,72380 

-  01 

■  9.2777D  -  01 

0.00000  -  00  - 

•6.71540  - 

01 

0.00000  -00 

O.OOOOO  -  00  ■ 

O.OOOOZ?  -  00 

9.81450  -01 

0.0000  0  - 

00 

0.00000  -  00 

0.00000  -  00 

1.3073Z)  -  02 

0.00000  -  00 

1.0229  0  - 

00 

0.00000  -  00 

O.OOOOO  -  00 

O.OOOOD  -  00 

21.44140  -  00 

0.0000  0  - 

00 

0.00000  -  00 

0.00000  -  00 

_  O.OOOOT?  -  00 

-6.28200  -  01 

O.OOOOO  - 

00 

0.00000  -  00 

0.00000  -  00 

8.4875Z)  -  04 

8.49750  -  04 

2.92290 

-03 

3.9602  0  - 

03 

1 

_ 

-1.4104Z)  -  05 

1.38490  -  05 

2.49550 

-  06 

-4.35480  - 

06 

9. 1794  0  -  nr, 

= 

3.1696D  -  05 

3.16960  -  05 

-1.27680 

-  05 

-2.14880  - 

05 

O.OOOOO  -  no 

2,5335D  -  05 

-2.40910  -  05 

9.58730 

-  05 

-1.64090  - 

04 

3.6092O  -  05 

1.2172D-03 

-1.19240  -  03 

1.46270 

-04 

-2.43660  - 

04 

1.36680  -  03 

■  -1.0088T9  -  03 

-1.00880  -  03 

-2.63490 

-03 

-3.54560  - 

03 

O.OOOOO  - 

00 

-1.3970D  -  05 

1.37180  -  05 

—  O' 

'.  jJ.'lO 

-  2 

'.98870  - 

05 

=  ' 

1.8899L>  -  05 

1.88990  -  05 

-5.04120 

-  05 

-7.23050  - 

05 

0,00000  - 

00 

-3.05200  -  04 

2.99690  -  04 

5.43020 

-  05 

-9.47500  - 

05 

-8.71390  - 

04 

8.94330  -  06 

-8.78180  -  06 

-1.59120 

-  06 

2.77650  - 

06 

2.5534  0  - 

05 
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Fable  B.3.  ACM  Entry  ARMA  Model  -  50%  Loss  of  Effectiveness  Left  Canard 

‘  -1.9488Z)  -  00  O.OOOOD-OO  0.3735D  -  01  O.OOOOD  -  00  UOnuOD  -  00 

O.OOOOZ)  -  00  -1.9845D-00  O.OOOOZ)  -  00  -101T5D-05  -1.0S805)-05 

.4i  =  -4.4830D  -  04  0,0000D-OO  -2.0113£)-00  D.OOOOD  -  00  0. 0000/9  -  00 

O.OOOOD  -  00  -24.84400-  00  0.00000  -  00  -9.77640  -  01  -6  J.53JO-01 

0.00000  -  00  5.75660  -  01  0.00000  -  00  -6.00590  -  04  -9,72000  -  01 

'  9.48850  -  01  0,00000  -  00  -4.96340  -  01  0,00000  00  O.OOOOO  -  nO  1 

0.00000  -  00  9.86660  -  01  0,00000  -  00  O.OOOOO  -  00  O.OOOOO  -  00 

.Ij  =  4,48200  -  03  0.00000  +  00  1. 00770  -  00  0,00000  -  uO  O.OOOOO-oO 

O.OOOOO  -  00  24.82600  -  00  O.OOOOO  -  00  0.00000  -  00  0.00000  -  00 

_  0,00000  -  00  -6,36240  -  01  0.00000  -  00  O.OOOOO  -  00  O.OOOOO  -  00  _ 

■-4,06290  -  03  -3,19700  -  03  5.98590  -  03  5,98590  -  03  O.OOOOO  -  00  ' 

-7.03210  -  06  -1.05790  -  05  3.99250  -  06  -3.99250  -  06  4,34370  -  05 

Ot  =  1.11910  -  05  2.23810  -  05  -2.17860  -  05  -2.17860  -  05  0.00000  -  00 

1.10820  -  05  -2.18090  -  05  1.64130  -  04  -1.64130  -  04  3,41960  -  05 

5.96190  -  04  -1.21800  -  03  1.94540  -  04  -1.94530  -  04  -1.16710  -  03  _ 

3  73230  -  03  2.78690  -  03  -5.17310  03  -5.4T34O-03  u.OOOOO  -  00  ' 

-7  0028O-  06  3.90100  -  05  3.97400  -  06  -3.97400  -  0»  -0.28910  -  06 

Ih  =  2,90580  -  05  3.60210  -  05  -4.80490  -  05  -4,80490  -  05  0  00000  -  00 

-1.76210  -  04  9  81580  -  04  9.99930  -  05  -9.99930  -  05  -158260  -  04 

4.51570  -  06  -2,51560  -  05  -2.56260  -  06  2.56260  -  6  1  05570  -  06  _ 
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Table  B.4.  ACM  Entry  ARMA  Model  -  25%  Loss  of  Effectiveness  Rudder 


-1.9684Z)  +  00 
0.0000/9  +  00 
-9,8466/9  -  03 
0,0000/9  4-  00 
0,0000/9  +  00 

9,68380  -  01 
0,0000D  +  00 
9,8443/9  -  03 
0,00000  +  00 
0,00000  +  00 


0,00000  +  00 
-1  99000  r  00 
0,00000  +  00 
-2,77790  +  01 
6,69700  -  01 

0,00000  +  00 
9,91890  -  01 
0,00000  +  00 
2,77500  +  0\ 
-7,19910  -  01 


4,35870  -  01 
0,00000  +  00 
-1,99240  +  00 
0,00000  +  00 
0,00000  +  00 

-4,10590  -  01 
0,00000  +  00 
9,84520  -  01 
0,00000  +  00 
0.00000  +  00 


0.00000  +  00 
-5,2378/9  -  06 
0.00000  +  00 
-9.75000  -  01 
-6.71510  -  04 

0.00000  +  00 
0.00000  +  00 
0.00000  +  00 
0.00000  +  00 
0.00000  +  00 


0  iinono  +  iin  1 
-1,2597/9  -  04  I 
li.m)[l0/9  -  iKi  ! 


-0,  u85i)/9  -  ill 
-9  -  11 ! 

i).nni)oo  +  oo 

O.OUUOO  +  00 
0,00000  4-  00 
D.OOOOO  -r  00 

o.ooooo  +  no 


7.31700  -  04 
8,28870  -  07 
3,17010  -  05 
2.33510  -  05 
-6.14140  -  07 

-7.45700  -  04 
8,26620  -  07 
2.41280  -  05 
2.31260  -  05 
-5.99960  -  07 


7.31700  -  04 
-8.28870  -  07 
3,17010  -  05 
-2.33510  -  05 
6  14140  -  07 

—  7.45700  —  04 
-8.26620  -  07 
2.41280  -  05 
-2.31260  -  05 
5.99960  -  07 


6.00400  -  03 
6.43520  -  04 
-2.18010  -  05 
1.38370  -  04 
-5  10260  -  02 

-3.81680  -  03 
6.43930  -  04 
-8.08470  -  05 
1.80150  -  02 
-4.67360  -  04 


6,00400  -  03 
-1.42150  -  05 
-2.18010  -  05 
-1.6374  0  -  04 
6.75650  -  04 

-5.81680  -  03 
-1.42040  -  05 
-8.08470  -  05 
-3.97390  -  04 
1.03100  -  05 


o.noooo  +  uo 

-1  2926  0  -  01 

0  ooili)//  III; 

!. 19700  -  i'5 
1  01170  -  ii2 

u.oouo/9  Oil 
-1,29350  -  04 
0.0000  0  +  on 
-3.61880  -  03 
9.38830  -  05 


j 
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Table  B.5.  ACM  Exit 

■-1.9837D-  00  -6.5693D-01 

2.6268Z)-07  -'.9658D-00 

.•ti  =  4.4540D-  04  1.4892D-04 

9.3315D-  08  -2.09520 -00 

-5.26660  -  09  1.17180  -  01 

9.83770  -  01  6.57250  -  01 

-2.60620  -  07  9.67960  -  01 

.42  =  _4.4191£)_  04  -1.48880  -  04 

-9.25830  -  08  2.09860  -  00 

5.22530  -  09  -1,15140  -Oi 

■  -2.43360  -  03  -2.45600  -  03 

1.54240  -  05  -1,89370  -  05 

fli  =  5.00310  -  06  4.84780  -  06 

5.66350  -  06  -1.23790  -  05 

3.32480  -  04  -3.02090  -  04 

2.25300  -  03  2.30310  -  03 

-1.89310  -  05  1.54620  -  05 

02  =  3.92680  -  06  3.75600  -  06 

-4.10490  -  05  3.35180  -  05 

2,25220  -  06  -1.83900  -  06 


ARMA  Model  -  No  Failures 

2.23860  -  00  1.16100  -  03  -186600  -  02  ' 

-1.18450  -  03  -2.53880  -  05  4  38270  -  04 

-1.97720  -  00  -1.3834  0  -  04  4,20160  -  06 

1.24400  -  03  -9.98310  -  01  -0,01100  -  02 

-6.65120  -  05  -9.18230  -  05  -9,98580  -  01  _ 

-2.23020  -  00  0,00000  -  00  0,00000  -  00  ' 

1.14760  -  03  0,00000  -  00  0.00000  -  90 

9.77580  -  01  O.OOOOO-OO  O.OOOOO  -  00 
1.23320  -  03  0.00000  -  00  O.OOOOD  -  00 

-6.80920  -  05  0.00000  -  00  0.00000  -  00  _ 

-2.33060  -  03  -2.33470  -  03  6.61960  -  07  ' 

2.C378O-  06  -2.03920  -u-3  1.75250  -  05 

-1.80590  -  06  -1.80590  -  06  1.58190  -  13 

1.96980  -  05  -1.96980  -  05  6.40600-06 

-3-35990  -  05  3,35990  -  05  -3.25660  -  04  _ 

2.42190  -  03  2,40330  -  03  -2.22570  -  06  ' 

1.66830  -  05  -1.66870  -  05  -2.27640  -  06 

-2-86210  -  06  -2.85880  -  06  6.55410  -10 

3.61740  -  05  -3,61760  -  05  -4  93540  -  06 

-1.98470  -  06  1.98480  -  06  2,70780  -  07 
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Table  B.6.  TF/TA  ARM  A  Model  -  No  Failures 


-.51.1370D  -  00 
O.OOOOD  -  00 
-1.8445D  -  01 
O.OOOOD  -  00 
O.OOOOD  -  00 

49.7810D  -  00 
O.OOOOD  -  00 
-7.6791D  -01 
O.OOOOD  -  00 
O.OOOOD  -  00 


O.OOOOD  -  00 
-'2.8733  D  -  05 
O.OOOOD  -  00 
-9.5265D  -  01 
-1.3'J54D  -  03 


O.OOOOD  -  00 
1.0815D  -  01 
O.OOOOD  -  no 
-  1.5014  D  -  00 
-9.2447D  -  01 


O.OOOOD  -  00 
O.OOOOD  -  00 
O.OOOOD  -  00 
O.OOOOD  -  00 
O.OOOOD  -  00 


O.OOOOD  -  no 
O.OOOOD  -  CO 
O.OOOOD  -  00 
O.OOOOD  -  00 
O.OOOOD  -  00 


-3.6847D  -  00 
O.OOOOD  -  00 
5.9207D  -  02 
O.OOOOD  -  00 
O.OOOOD  -  00 

2.6832D  -  00 
O.OOOOD  -  00 
-5.9154D  -  02 
O.OOOOD  -  00 
O.OOOOD  -  00 

1.0618D  -  02 
-1.8033D  -  05 
7.3451D  -  05 
7.1020D  -  05 
3.9102D  -  03 

2.n494D  -  02 
-7.5883D  -  05 
-4.2691D  -  04 
-4.6779D  -  03 
1.2979D  -  04 


O.OOOOD  -  00 
-1.9519D  -  00 
O.OOOOD  -  00 
-59.2300D  -  00 
1.4404D  -  00 

O.OOOOD  -  00 
9.5785D  -  01 
O.OOOOD  -  00 
59. 0480 D  -00 
-1.6383D  -  00 

-1.0618D  -  02 
1.8033D  -05 
7.3451D  -  05 
-7.1020D  -  05 
-3.9in2D  -  03 

2.0494D  -  02 
7.5883D  -  05 
-4.2691D  -  04 
4.6779D  -  03 
-1.2979D  -  04 


1.5259D  -  02 
-4.1580D  -  06 
-4.0607D  -  05 
2.9993D  -  04 
1.0870D  -  04 

-3.6496  D  -  02 
1.0468D  -  05 
7.9098D  -  04 
6.4531  D  -  04 
-1.7904D  -  05 


1.5259D  -  02 
9.1373D  -  06 
-4.0607D  -  05 
-2.9997D  -  04 
-1.0819D  -  04 

-3.6496D  -  02 
-1.5268D  -  05 
7.9098D  -  04 
-9.4120D  -  04 
2.6114D  -  05 


O.OOOOD  -  00  ■ 
1.0194D  -  04 
O.OOOOD  -  00 
1.0971D  -  04 
-4.4117D  -  03 

O.OOOOD -no 
1.2956D  -  05 
O.OOOOD  -  00 
7.0867D  -  04 
-2.2160D  -  05 


Table  B.7.  TF/TA  ARM  A  Model  -  30%  Loss  of  Effectiveness  Left  Trailing  Edge 

■-1.9289D-00  O.OOOOD-  00  2.4122Z?  -  00  O.OOOOD  -  00  O.OOOOD  -  00  ‘ 

O.OOOOD-  00  -1.956l£)-  00  0.0000£)  -  00  -2.5617D  -  05  1.3755D  -  05 

Ai  =  7.237lZ)-03  O.OOOOD  -  00  -1.9644£)  -  00  O.OOOOZ9  +  00  O.OOOOD  -  00 

O.OOOOD-  00  57.5160£)-00  O.OOOOZ)  -  00  -9.5402D  -  01  -1.4610£)  -  00 

O.OOOOD -  00  1.6261D-00  O.OOOOD  -  00  -1.5406D-  03  -9.2005D  -  01 

9.2900D-01  O.OOOOD -  00  -2.3558D  -  00  O.OOOOD  -  00  0.0000D-00‘ 

O.OOOOD -  00  9.6193D-01  O.OOOOD  -  00  O.OOOOD  -  00  O.OOOOD  -  00 

.-Ij  =  -7.2311D-03  O.OOOOD -  00  9.7021D-01  O.OOOOD  -  00  O.OOOOD  -  00 

O.OOOOD -  00  57.7350D-00  O.OOOOD  -  00  O.OOOOD  -  00  O.OOOOD -00 

O.OOOOD -  00  -1.8209D-00  O.OOOOD -00  O.OOOOD  -  00  O.OOOOD  -  00  , 

■-8.8814D-  03  -8.8814D-  03  1.0201D-  02  1.3521D  -  02  O.OOOOD  -  00  ‘ 

-2.2384D-0S  2.2384D  -  05  2.2204D  -  06  -3.6512D-  06  1.0626D  -  04 

=  7.2684D-  05  7.2684D-  05  -2.3547D  -  05  -3.9662D  -  05  O.OOOOD  -  00 

7.0004D-  05  -7.0004D-  05  1.7183D  -  04  -2.9585D  -  04  1.1159D-04 

3.8607D-  03  -3.8607D  -  03  4.3589D  -  05  -8.8733D  -  05  -4.3508D  -  03 

■  8.2180D-  03  8.2180D-  03  -9.3663D  -  03  -1.2370D-02  O.OOOOD  -  00  ‘ 

-7.0757D-  05  7.0757D-  05  2.1920D  -  06  -3.6047D  -  06  7.5178D  -  06 

=  7.4294D  -  06  7.4294D-  06  5.0055D  -  05  5.7807D  -  05  O.OOOOD  -  00 

-4.2185D-03  4.2185D  -  03  1.3069D  -  04  -2.1491D-  04  4.4821D-U4 

1.3394D-  04  -1.3394D-  04  -4.1495D  -  06  6.8238D  -  06  -1.1231D-05 
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Table  B.8.  TF/TA  ARMA  Model  -  50%  Loss  of  Effectiveness  Left  Canard 

■  -1.9180D  -  00  0.0000£)-  00  2.4247Z)  -  00  O.OOOOZ)  -  00  O.OOOOZ)  -  00 

O.OOOOD  -  00  -1.9713D-  00  0.0000 £>-  00  -1.2614Z)-  05  -5.2084D-04 

.-1,  =  4.7557D-03  O.OOOOZ) -00  -1.9509Z)-00  O.OOOOZ)  -  00  O.OOOOZ) -00 

O.OOOOZ) -  00  -83.6540D-00  O.OOOOZ)  -  00  -9.3315Z)  -  01  -2.12890-00 

O.OOOOD  -  00  2.0604Z)-00  O.OOOOD  -  00  -1.9038Z)-  03  -9.033nZ)  -  01 

9.i806Z)-oi  O.OOOOZ) -  00  -2.3596Z)-oo  O.OOOOZ)  -  00  o.oonnZ)-no' 

O.OOOOZ) -  00  9.7736Z)-01  O.OOOOZ) -00  O.OOOOZ)  -  DU  O.onoOZ)  -  Oil 

.li  =  -4.7512Z)-03  O.OOOOZ) -  00  9.5472Z)-01  O.OOOOZ)  -  00  O.OOOOZ)  -  DO 

O.OOOOZ) -  00  83.3500D-00  O.OOOOZ) -00  O.OOOOZ)  -  00  O.OOOOD  -  00 

O.OOOOD -  00  -2.2661D-00  O.OOOOD -00  O.OOOOD  -  00  O.OOOOD  -  00  _ 

■  -1.3180D-  02  -9.7372D-03  1.6645D-  02  1.6645D-02  O.OOOOD -00' 

-3.0171D-  06  5.2788D-  06  3.0283D  -  06  -3.0283D  -  06  1.1496D-04 

Di  =  4.2701D-  05  8.4482D-05  -4.6574D-  05  -4.6574D-05  O.OOOOD  -  00 

3.6503D-  05  -8.5166D-  05  3.4423D  -  04  -3.4423D  -  04  1.2590D-04 

2.2580D-  03  -4.4667D-  03  6.4002D  -  05  -6.4002D  -  05  -5  0719D  -  03  _ 

1.2025D-  02  8.7762D-  03  -1.4993D-  02  -1.4993D-02  O.ooiidZ)  -  DO  ' 

-5.2586D-  05  1.0441D-  04  2.9923D  -  06  -2.9923D  -  06  l.O.DOiiD  -  0-5 

D2  -  -2.0516D-  05  3.7108D-  05  3.2466D  -  05  3.2466D  -  05  O.OOOOD  -  00 

-4.4845D-  03  8.9044D  -  03  2.5519D  -  04  -2.5519D  -  04  l.  tl53D-03 

1.2193D-  04  -2.4210D-  04  -6.9381D-  06  6.9381D-  06  -3  8l.8lD-05_ 


B-19 


Table  B.9.  TF/TA  ARMA  Model  -  25%  Loss  of  Effectiveness  Rudder 


-1.91'29Z)  +  00 

0.00000  +  00 

1.71080  +  00 

0.00000  +  00 

0.00000  +  00 

O.OOOOD  +  00 

-1.98280  +  00 

0.00000  +  00 

-5.69970  -  06 

-5.66670  - 

04 

.-ll  = 

-1.2382Z)  -  02 

0.00000  +  00 

-1.95450  +  00 

0.00000  +  00 

0.00000  +  00 

O.OQQOD  +  00 

-6.66910  +  01 

0.00000  +  00 

-9.46700-  01 

-1.69080  +  00 

O.OOOOD  +  00 

2.33550  +  00 

0.00000  +  00 

-2.04710  -  03 

-9.05180  - 

-01  J 

■  9.12980-01 

0.00000  +  00 

-1.64140  +  00 

0.00000  +  00  0.00000  +  00  ■ 

0.00000  +  00 

9.87170  -  01 

0.00000  +  00 

0.00000  +00  0.00000  +  00 

'2  = 

1.23710  -  02 

0.00000  +  00 

9.44540  -  01 

0.00000  +  00  0.00000  +  00 

0.00000  +  00 

6.64750  +  01 

0.00000  +  00 

0.00000  +  00  0.00000  +  00 

_  0.00000  +  00 

-2.47770  +  00 

0.00000  +  00 

0.00000  +  00  0.00000  +  00  _ 

-8.36080  -  03 

-8.36080  -  03 

1.39060  -  02 

1.39060  -  02 

O.OOOOO  +  00 

2.59120  -  05 

-2.59120  -  05 

1.89440  -  03 

-1.63500  -  05 

-3.49260  - 

04 

!h  = 

7.22960  -  05 

7.2296  0  -  05 

-3.91500  -  05 

-3.91500  -  05 

O.OOOOO  +  00 

6.86010  -  05 

-6.86010  -  05 

2.81360  -  04 

-2.95900  -  04 

8.52610  - 

05 

-7.07560  -  07 

7.07560  -  07 

-1.50690  -  01 

9.70930  -  04 

2.99220  - 

02 

7.53510  -03 

7.53510  -  03 

-1.25520  -  02 

-1.25520  -  02 

0.00000  +  00  ■ 

-2.37350  -  05 

2.37350  -  05 

1.91520  -  03 

-1.64630  -  05 

-4.0302O  - 

04 

= 

1.74080  -04 

1.7408  0  -  04 

-2.0952O  -  04 

-2.09520  -  04 

O.OOOOO  +  00 

-1.59830  -  03 

1.59830  -03 

1.28970  -  01 

-1.10860  -  03 

-2.71390  - 

02 

5.95720  -  05 

-5.95720  -  05 

-4.80680  -  03 

4.13210  -  05 

1. 01150  - 

03 
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Appendix  C.  GAIN  MATRICES  AND  STEP  RESPONSE  ROOTS 


This  Appendix  contains  the  fixed-gain  matrices  for  the  no  failure  flight  conditions 
and  the  wi  roots  for  the  step  response  PI  controller  design. 
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Table  C.l.  Discrete  PI  Controller  Gain  Matrices  -  ACMENTRY 


■  0.2750 

2.6579 

22.6167 

0.9347 

3.2537  ■ 

0.2750 

-2.6579 

22.6167 

-0.9347 

-3.2537 

= 

0.3989 

-0.8640 

-2.9566 

2.4292 

-0.1156 

0.3939 

0.8640 

-2.9566 

-2.4292 

0. 1 156 

_  0.0000 

4.5235 

0.0000 

2.3043 

-3.4078  _ 

■  0.0825 

0.7974 

6.7850 

0.2804 

0.9761  ■ 

0.0825 

-0.7974 

6.7850 

-0.2804 

-0.9761 

«2  = 

0.1197 

-0.2592 

-0.8870 

0.7288 

-0.0347 

0.1197 

0.2592 

-0.8870 

-0.7288 

0.0347 

0.0000 

1.3570 

0.0000 

0.6913 

-1.0223 

Table  C. 

2.  Discrete  PI  Controller  Gain  Matrices  -  ACMEXIT 

■  -0.2877 

1.8139 

118.6223 

5.9943  7.1919  ‘ 

-0.2728 

-1.8515 

112.4901 

-6.1188  -7.3413 

A'l  = 

-0.8052 

-1.4469 

-99.0041 

19.5178  0.7930 

-0.7271 

1.5019 

-131.2262 

—  19.3359  —0.5747 

_  -0.0326 

3.8739 

13.4352 

7.9424  -24.40' 

■4  J 

■  -0.0863 

0.5442 

35.5867 

1.7983  2.1576  ' 

-0.0818 

-0.5555 

33.7470 

-1.8356  -2.2024 

I<2  = 

-0.2416 

-0.4341 

-29.7012 

5.8553  0.2379 

-0.2181 

0.4506 

-39.3679 

-5.8008  -0.1724 

-0.0098 

1.1622 

4.0306 

2.3827  -7.3222 
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Table  C.3.  Discrete  PI  Controller  Gain  Matrices  -  TF/TA 


A'l 


0.1372 

0.8360 

15.3257 

0.1741 

0.8983 

0.1549 

-0.8360 

15.9544 

-0.1741 

-0.8983 

= 

0.2707 

-0.4656 

9.5751 

1.3828 

0.0093 

0.2608 

0.4656 

9.2249 

-1.3828 

-0.0093 

_  -0.0155 

1.4562 

-0.5475 

0.3824 

-1.1936 

0.0411 

0.2508 

4.-5977 

0.0522 

0.2695 

0.0465 

-0.2508 

4.7863 

-0.0522 

-0.2695 

0.0812 

-0.1397 

2.8725 

0.4148 

0.0028 

0.0782 

0.1397 

2.7675 

-0.4148 

-0.0028 

-0.0046 

0.4368 

-0.1643 

0.1147 

-0,3581 

Table  C.4.  Step  Response  PI  Controller  Gain  Matrices  -  ACM  Entry 


5.4705 

319.1986 

145.8124 

13.6863 

39.6138 

5.4705 

-319.1986 

145.8124 

-13.6863 

-39.6138 

7.9174 

-103.6666 

-18.4503 

143.8869 

-7.2031 

7.9174 

103.6666 

-18.4503 

-143.8869 

7.2031 

0.0000 

-543.6125 

0.0000 

61.1122 

13.4587 

0.0028 

8.4049 

6.9320 

-3.8274 

106.9155 

0.0028 

-8.4049 

6.9320 

3.8274 

-106.9155 

0.0027 

1.3925 

7.9268 

1.2441 

-34.6890 

0.0027 

-1.3925 

7.9268 

-1.2441 

34.6890 

0.0000 

24.4185 

0.0000 

-6.5138 

180.3718 
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Table  C.5.  Step  Response  FI  Controller  Gain  Matrices  -  ACM  Exit 


-0.0058 

0.2178 

0.7526  - 

■0.0085 

0.0583  ■ 

-0,0055 

-0.2223 

0.7138 

0.0093  - 

•0,0596 

-0.0163 

-0.1734 

-0.6660 

1.3210  - 

■0.0077 

-0.0147 

0.1801 

-0.8705  - 

■1.3197 

0.0095 

-0.0007 

0.4651 

0.085*  - 

■0.2897  - 

■0.1137 

-0.8029 

5.4103 

-94.2627 

9.1723 

85.0569 

-0.7614 

-4.9904 

-82.6466 

13.2984 

—  78.4553 

-1.6540 

9.2572 

-183.2629 

27.4612 

-59.1481 

-1.4359 

-8.1117 

-164.1554 

20.5369 

77.1570 

-0.0909 

7.0755 

-17.6024 

-3.6862 

174. 80i 
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Table  C.6.  Step  Response  PI  Controller  Gain  Matrices  -  TF/TA 


2.7661 

100.5031 

108,0596 

4.5846 

11.7829 

3.1223 

-100.5031 

113,2086 

-4.5846 

-11.7829 

5.4248 

-55.8763 

78.4220 

80.0561 

3.4695 

-0.3102 

175.4099 

-4.4842 

12.0628 

3.7161 

0.0411 

0.2508 

4.5977 

0,0522 

0.2695  ■ 

0.0465 

-0.2508 

4.7863  - 

-0.0522  - 

•0.2695 

0.0812 

-0.1397 

2.8725 

0.4148 

0.0028 

0.0782 

0.1397 

2.7675  - 

•0.4148  - 

■0,0028 

-0.0046 

0.4368  - 

■0.1643 

0.1147  - 

•0.3581 
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Table  C. 


7.  w/  plane  roots  for  ACM  Entry  and  TFTA  for  the  step  response  method 


ACM  ENTRY 

TFTA 

A.TOE-'^ 

-4.70£-'^ 

6.G3E-® 

6.60f;-*^ 

-2.34E-* 

-2.45E-* 

-2.86E-3±;2.10E-3 

-3.73£-3  ±;3.18£-^ 

-2.2034  ±j  .7398 

-2.3131 

-2.8844  ±j  1.3522 

-2.5416 

-1.4674  ±j  4.8738 

-8.1269 

-9.9482  (5) 

-2.8595  ±j  8.5404 

-1.1482  ±;  11.652 

-9.9482  (5) 

-45.617  ±j  27.997 

-11.434 

-55.799  ±]  31.105 

-16.310 

-56.411  ±]  31.192 

-24.174  ±;  44.036 

-56.522  ±;  31.006 

-51.822 

-57.916  ±j  31.327 

-50.148  ±j  30.743 
-55.239  ±j  29.878 
-54.500  ±]  31.901 
-.57.237  ±]  31.264 
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Appendix  D.  ADAPTIVE  ALGORITHM 


The  following  source  code  for  the  adaptive  algorithm  is  a  FORTRAN  subroutine 
obtained  from  Hammond  (4)  which  has  to  be  linked  to  MATRIX\.  This  algorithm 
will  require  more  utfaui..  memory  than  found  in  M AT RI X x ^  so  the  st.ick  size  must  be 
increased  in  the  MATRIXx  FORTRAN  code.  Specific  instructions  to  increase  the  stack 
size  and  link  FORTR.A..N'  code  can  be  found  in  the  M ATRIXx  users  manual.  .Also  notc 
that  the  most  recent  versions  of  the  IMSL  Library  and  M AT RI X x  are  used. 


C - - - 

Cl  I 

Cl  THIS  PROGRAM  PROVIDES  A  RECURSIVE  LEAST  SQUARES  ESTIMATION  I 
Cl  OF  THE  PARAMETERS  FOR  THE  CRCA  5  X  5  A  AND  B  COEFFICIENTS.  I 
Cl  LAST  REVISION  11  SEPTEMBER  1989  I 

Cl  THIS  PROGRAM  SHOULD  BE  LINKED  TO  VERSION  10.0  IMSL  LIBRARY  I 
Cl  AND  WAS  WRITTEN  TO  RUN  ON  MATRIXX  VERSION  7.0.  I 

Cl  SEE  till  FOR. COM  FOR  THE  APPROPRIATE  COMMANDS  I 

C . . 

SUBROUTINE  UPDUSR  (INFO .NUMBER, T.U.NU, X ,XD0T, NX ,Y , NY , 

+  RP.IP) 

DOUBLE  PRECISION  T,U(*) ,X(*) ,Y(*) ,XD0T(*) ,RP(*) 

INTEGER  INF0(4),  IP(*) 

CHARACTER*3  CNUM 
C 

INTEGER  MAXNUM 
PARAMETER  (MAXNUM=1) 

C 

IF  ( NUMBER. GT. MAXNUM)  THEM 
INF0(l)*-2 

WRITE(CNUM,111)  NUMBER 
111  F0RMAT(I3) 

CALL  MATWRC  ’) 

CALL  MATWR 

+  (’SIM.ERROR:  NOT  ABLE  TO  UPDATE  USER  FUNCTION ’ //CNUM) 

RETURN 
ENDIF 

CALL  USROl  (INFO,T,U,NU,X,XDOT,NX,Y,NY.RP,IP) 

RETURN 

END 
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c 

c 

c 

c 

c 


!  START  OF  USER  SUBROUTINE  FOR  THE  RECURSIVE  LEAST  SQUARES 
I  ESTIMATION  AND  GAIN  MATRIX  CALCULATIONS  FOR  THE  ARMA 
I  MODEL  REPRESENTATION 


SUBROUTINE  USROl (INFO ,T .U ,NU.X .XDOT.NX , Y , NY ,RP , IP) 
DOUBLE  PRECISION  T.U(+) .X(*) ,XDOT(*) , Y(*) .RP(*) , 


+ 


+ 


INTEGER 

LOGICAL 


XNPLUSdOO.S)  ,GAMD1(5,5)  .XNPLUST(5 , 100) , 
XNPLUSPN(5,100) ,PN(100,100) ,GAMMA(5,5)  , 
F0RGET(5,S) ,PNXNPLUS(100,5) , 

G0N(5,5)  ,HINV(5,5)  ,A1ARMA(5 ,5pi  i(2)  . 
A2ARMA(5,5) ,B1ARMA(5,5) ,B2ARMA(5,5) , 
KKS.S)  ,K2(5,5)  .PND2(100,100) , 

THED2(5) .THED4(100) , 

THETA(IOO) ,EYE5(5,5) ,GAMMAI(5,5) ,THED3(5) , 
SIGMA(5) ,G0(5.5) ,G0DEN(5,5) ,PIT(5) , 
G0DENI(5,5) .PNPLUSClOO.lOO) ,THETA1(100)  , 
PXPLGAMClCO.S) .WK4REA(40) ,Gni(5,5) 

IP(*)  .  INF0(4) ,N,IA,IDGT,IER,I,J,K.L 
INIT, STATE, OUTPUT 


INIT  =  INFO(2).NE.O 

STATE  =  INFO (3) .NE.O 

OUTPUT^  INFO (4). NE.O 

IF  (STATE  .OR.  (.NOT.  OUTPUT))  THEN 

GOTO  999 

Ftjnrr 

N=5  ! PARAMETERS  FOR  IMSL  LIBRARY 

LDA=5  ! PARAMETERS  FOR  IMSL  LIBRARY 

LDAINV=5  ! PARAMETERS  FOR  IMSL  LIBRARY 

C - - - - - - - . - 

Cl  CHECK  FOR  INITIALIZATION  OF  THE  THETA(O)  AND  P(0) 

Cl - - - - - - - 

TSTART».050  !3  PERIODS 

IF  (TSTART.GT.T)  THEN 

OPEN (UNIT= 102 , FILE= ’ THETANOM . RLS ’ . STATUS= ' OLD ’ ) ! THETA 
DO  1=1,100 

READ ( 102, =)THETA(I) 

END  DO 
CLOSE(102) 

OPEN(UNIT=103,FILE=*PONOM.RLS’ ,STATUS=’OLD’ )  !PN 
DO  1=1,100 
DO  J=l,100 

READ(103,*)PN(I,J) 

END  DO 
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END  DO 
CL0SE(103) 

GOTO  600  ! INSURES  FULLY  POPULATED  XNPLU5 
ENDIF 


C 

n 

c 


CALCULATION  OF  NEW  PARAMETER  VECTOR 


CALL  MMULTl 
DO  1=1,5 
THED3(I) 
END  DO 
CALL  MMULTl 
DO  J=l,100 
THETAl(J) 
END  DO 


(5,100,1, XNPLUST , THETA . THED2 ) 
=  U(I+10)-THED2(I) 

(100,5,1, PXPLGAM ,THED3 ,THED4) 
=  THETA(J)  ♦  THED4(J) 


C . . - - - 

Cl  UPDATE  PARAMETER  VECTOR  FOR  NEXT  ITERATION 
C . - . - . -  — . - 

DO  :.=  i,ioo 

THETA(I)=THETA1(I) 

END  DO 


C . - . 

Cl  '"HE  RECURSIVE  UPDATE  EQUATION  STARTS  HERE.  I 

C . 

C . - . 

Cl  UPDATE  XNPLUS  WITH  THE  OUTPUT  FOR  NEXT  ITERATION  I 

C . . . . . - . 

600  K=1  ! COUNTER  FOR  UPDATE  MATRIX 

DO  J=75,95,5 
DO  1=1,5 

XNPLUS(J+I,K)  =  U(I+5)  !U(I+5)--U(T-2) 

END  DO 
K=K  +  1 
END  DO 


K=1  ! COUNTER  FDR  UPDATE  MATRIX 

DO  J=50,70,5 
DO  I  =  1,5 

XNPLUS(J+I,K)=U(I)  !U(I)— U(T-l) 

END  DO 
K=K  +  1 
END  DO 


C 


K=1 

DO  J=25,45,5 
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DO  1=1,5 

XNPLUS(J+I ,K)  =  -0(1+15)  lUPDATES  Y(T-2) 
END  DO 
K=K+1 
END  DO 
K=1 

C . . - - - - -  - 

DO  J=0,20,5 
DO  1=1,5 

XNPLUS(J+I ,K)=  -U(I+10)  !  UPDATES  Y(T-l) 

END  DO 
K=K+1 
END  DO 


Cl  CALCULATION  OF  P(K)  ♦  XN(K+  1) 


CALL  MMULT  ( 100 , 100 , 5 ,PN . XNPLUS .PNXNPLUS) 

r 

Cl  CALCULATE  XNPLUST  *  PN 

C . . . . . . . 

C- . - . - . - . 

Cl  THE  FOLLOWING  LOOP  UPDATES  THE  FORGETTING  FACTOR  MATRIX 

C . - . 

DO  1=1,5 
DO  J=l,5 

IF(I.Eq.J)THEN 

F0RGET(I,J)=U(31) 

ELSE 

FORGETCI, J)=0 
END  IF 
END  DO 
END  DO 
DO  1=1,5 

DO  J=l,100 

XNPLUST ( I, J)  =  XNPLUS (J, I) 

END  DO 
END  DO 

CALL  MMULT  (5 , 100 , 100 .XNPLUST, PN .XNPLUSPN) 

C--- . - - - - - 

Cl  CALCULATE  X(N+1)T*PN*X(N+1) 


CALL  MMULT  (5 , 100 , 5 , XNPLUSPN , XNPLUS .GAMDl ) 


Cl  CALCULATE  GAMMA (N+1) 

C- . - - - 
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DO  1=1,5 
DO  J=1.5 

GAMMAd.J)  =  FORGETCI.J)  +  GAMDl(I.J) 

END  DO 
END  DO 

C--- - - - - - - - - - 

Cl  COMPUTE  GAMMA (N+1)  INVERSE  USING  IMSL  ROUTINE 

C - - - - - - - . — . . 

CALL  DLINRG  (N .GAMMA , LDA .GAMMAI .LDAINV) 

CALL  MMULT  ( 100 , 5 ,5 .PNXNPLUS .GAMMAI .PXPLGAM) 

C - - - 

Cl  CALCULATE  PNPLUS.  UPDATE  OF  COVARIANCE  MATRIX 

C . . . . . . . 

CALL  MMULT  (100 .5 . 100 .PXPLGAM.XNPLUSPN.PND2) 

C  =  =  =  =  ===s==s=  =  =  =  =  r  =  =  =:  =  =  =  ===  =  =  =  =  =  »==.=====  =  =  ===  =  =  s  =  =  =  =  =  =  =  =  =  =  = 

DO  J=1.100 
DO  K=l,100 

PNPLUS(J.K)=PN(J,K)-PND2(J.K) 

END  DO 
END  DO 

C . . . 

Cl  UPDATE  COVARIANCE  MATRIX  AND  PREPARE  FOR  NEXT  ITERATION  I 

Cl- . - . 

DO  1=1,100 
DO  J=l,100 

PN(I,J)  =  PNPLUSCI.J) 

END  DO 
END  DO 

C - - - - - 

Cl 

Cl  THIS  CONCLUDES  THE  BASIC  CALCULATION  OF  THE  LEAST  SQUARES 
Cl  ESTIMATE  OF  THE  PARAMETER  VECTOR  THETA.  NEXT  THE  GAIN 
Cl  MATRICES  WILL  BE  CALCULATED 

Cl 

C - - - . . . . . - . - 

C . - . - - - 

Cl 

Cl  CALCULATION  OF  THE  IDENTITY  MATRIX  FOR  5X5 

C-- . - - - - - - 

DO  1=1,5 
DO  J=l,5 

IFCI.EQ. J)THEN 
EYE5(I,J)=1 
ELSE 

EYE5(I,J)=0 
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END  IF 
END  DO 
END  DO 

C - CALCULATION  OF  H  AND  H  INVERSE . - . . . . 

K=0 

DO  1=1,5 
DO  J=l,5 

A1ARMA(I,J)=THETA(K  +  J) 

A2ARMA(I, J)=THETA(2S  +K  +  J) 

B 1 ARMA ( I , J ) =THETA ( 50+K+ J) 

B2ARMA(I.J)=THETA(75+K+J) 

GON(I.J)  =B1ARMA(I.J)  +  B2ARMA(I.J) 

GODEN(I,J)  =EYE5(I.J)  +  A1ARMA(I,J)  +  A2ARMA(I,J) 
END  DO 
K=K  +  5 
END  DO 

CALL  DLINRG  (N. BURMA, LDA.HINV.LDAINV) 

CALL  DLINRG  (M,GODEN,LDA,GODENI,LDAINV) 

C . - . 

Cl  INITIALIZE  SIGMA  MATRIX 

C--- . 

DO  1=1,5 

SIGMA(I)=U(I+20) 

END  DO 

C . 

Cl  INITIALIZE  PI  MATRIX 

C . . - . . 

DO  1=1,5 

PIT(I)=U(25+I) 

END  DO 

C - - - - - - - 

Cl  CALCULATION  OF  G(0)  AND  INV(GO) 

C--- . - - - - - 

CALL  MMULT  (5 ,5 ,5 ,GODENI ,GON,GO) 

CALL  DLINRG  (N ,GO ,LDA ,GOI ,LDAINV) 

C . . — . . 

Cl  CALCULATION  OF  K1  AND  K2  FOR  GAIN 

C . - . - . 

DO  1=1,5 
DO  J=l,5 

K1(I,J)=  HINV(I,J)*SIGMA(J) 

K2(I,J)=  GOI(I,J)*PIT(J) 

END  DO 
END  DO 
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c - - - - - - 

Cl  THE  FOLLOWING  ROUTINES  WILL  PREPARE  THE  OUTPUTS  OF  THE  USER  I 
Cl  BLOCK.  Y(l)  -  Y(25)=K1,  Y(26)-Y(50)=K2,Y(51)-Y(150)=THETA  I 

Cl- . - . - . . . . 

L=25 

K=0 

DO  1=1,5 
DO  J=l,5 

Y(J+K)=K1(I,J) 

Y(J>L)=K2(I,J) 

END  DO 
K=K+5 
L=L+5 
END  DO 
DO  1=1,100 

Y(I+50)=THETA(I) 

END  DO 
999  RETURN 
END 

. . 

Cl  SUBROUTINE  MMULT  FOR  THE  MATRIX  MULTIPY-SIZE  >  1 

C . 

SUBROUTINE  MMULT  (MATARW,MATACL,MATABCL,MATA,MATB,MATC) 
DOUBLE  PRECISION  MATC(MATARW,MATABCL) ,MATA(MATARW ,MATACL) , 
+  MATB(MATACL,MATABCL) 

INTEGER  MATARW,MATACL,MATABCL 
DO  I=1,MATARW 
DO  K=1,MAT.'IBCL 

SUM=0  ! RESET  PRODUCT  SUM 

DO  J=1,MATACL 

SUM=SUM  +  MATA(I,J)*MATB(J,K) 

END  DO 

MATC(I,K)=SUM 
END  DO 
END  DO 
RETURN 
END 

C . . . . . 

Cl  MATRIX  MULTIPLY  ROUTINE  FOR  MATRICES  WITH  A  COLUMN 

C . . . . . . 

SUBROUTINE  MMULTl  (MATARW,MATACL.MATABCL,MATA,MATB,MATC) 
DOUBLE  PRECISION  MATC(MATARW) ,MATA(MATARW,MATACL) , 

+  MATB(MATACL) 

INTEGER  MATARW,MATACL,MATABCL 
DO  I=1,MATARW 
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! RESET  PRODUCT  SUM 


SUM=0 

DO  J=1,MATACL 

SUM  =  SUM  +  MATA(I,J)*MATB(J) 
END  DO 

MATC(I)  =  SUM 
END  DO 
RETURN 
END 
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